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ABSTRACT 

Instabilities driven by thermal and lepton diffusion (doubly diffusive instabilities) in a Ledoux stable 
fluid will, if present below the neutrinosphere of the collapsed core of a supernova progenitor (proto- 
supernova), induce convective-like fluid motions there. These fluid motions may enhance the neutrino 
emission by advecting neutrinos outward toward the neutrinosphere, and may thus play an important 
role in the supernova mechanism. "Neutron fingers," in particular, have been suggested as being crit- 
ical for producing explosions in the soph isticated spherically symmetric supernova simulations by the 
Liver more group (I Wilson fc Mavlelll993l e.g.,). These have been argued to arise in an extensive re- 
gion below the neutrinosphere of a proto-supernova where entropy and lepton gradients are stabilizing 
and destabilizing, respectively, if, as they assert, the rate of neutrino- mediated thermal equilibration 
greatly exceeds that of neutrino-mediated lepton equilibration. Application of the Livermore group's 
criteria to models derived from core collapse simulations using both their equation of state and the 
Lattimer-Swesty equation of state do indeed show a large region below the neutrinosphere unstable 
to neutron fingers. 

Because of the potential importance of fluid instabilities for the supernova mechanism, and the 
desire to understand the origin of convective-like fluid motions that may arise in upcoming multi- 
dimensio nal radiation- hydro dynam ical simulations of core collapse, we develop a methodology intro- 
duced bv iBruenn Dineval (|1996T ) for analyzing the stability of a fluid in the presence of neutrinos 
of all flavors and in the presence of a gravitational field. Neutrino-mediated thermal and lepton equi- 
libration between a fluid element and its surroundings (background) is modeled as a linear system 
characterized by four response functions (i.e., thermal and lepton equilibration driven by entropy and 
lepton fraction differences between a fluid element and the background), the latter evaluated for a 
given thermodynamic state and fluid element radius by detailed neutrino transport simulations. These 
transport simulations employ both traditional and improved neutrino physics. 

When applied to an extensive two-dimensional grid of core radii and fluid element sizes for each 
of several time slices of a number of proto-supernovae, we find no evidence for the neutron finger 
instability as described by the Livermore group. We find, instead, that the rate of lepton equilibration 
always exceeds that of thermal equilibration. Furthermore, we find that neither of the "cross" response 
functions, that is, entropy equilibration driven by a lepton fraction difference, and lepton equilibration 
driven by an entropy difference, is zero and that the first of these tends to be the largest of the four 
response functions in magnitude. These cross response functions play a critical role in the dynamics 
of the equilibration of a fluid element with the background. An important consequence of this is 
the presence of a doubly diffusive instability, which we refer to as "lepto-entropy fingers," in an 
extensive region below the neutrinosphere where the lepton number, Y^, is small. This instability 
is driven by a mechanism very different from that giving rise to neutron fingers, and may play an 
important role in enhancing the neutrino emission. Deep in the core where the entropy is low and the 
lepton number higher, our analysis indicates a region unstable to another instability, also involving 
the cross response functions, which we refer to as "lepto-entropy semiconvection." These instabilities, 
particularly lepto-entropy fingers, may have already been seen in some multi-dimensional core collapse 
simulations described in the literature. 

Subject headings: (stars:) supernovae: general - neutrinos - fluid instabilities 



1. INTRODUCTION 

It has long been recognized that fluid instabilities 
may play an important role in the supernova mech- 
anism. It is generally agreed that entropy driven 
convection between the neutrinosphere and the shock 
during the postbounce phase of the collapsed core of 
a supernova progenitor (hereafter "proto-supernova") 



will contribute turbulent pressure and increa s e the 
neutrino en ergy-deposition effi c iency (iBethel [1 990; 
MilleTTnn ITflB iHepnt etall ITflfll iRnrrows et »,1J 
1995; Janka fc Mullerl 11996c i Mezzacapoa et alJ Il998b: 
lFrverl1999UFrver fc Warreul2002HBuras et all200fl. but 
this may not be enough to produce explosions. 

More controversial is the nature and the role of fluid 
instabilities below the neutrinosphere. Here matter and 
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neutrinos are tightly coupled, and convective-like mo- 
tions induced by instabilities can advect neutrinos. If 
these convective-like motions occur between the deeper 
interior and the neutrinosphere, i/qS could be advected 
from the lepton rich deeper interior to the neutri- 
nosphere, thus adding to the rate of v e diffusion and 
enhancing the u e luminosity. The analysis of fluid in- 
stabilities below the neutrinosphere must take into ac- 
count the possibilities of thermal and lepton transfer by 
neutrinos. For example, a displaced fluid element in the 
presence of an entropy and/or lepton fraction gradient 
will find that its entropy and/or lepton fraction will dif- 
fer from that of its surroundings or background. This 
will induce neutrino mediated thermal and lepton trans- 
fer which can modify the hydrodynamic buoyancy forces 
felt by the fluid element. This in turn can modify the 
subsequent growth rate of an instability, or the very ex- 
istence of an instability. In particular, instabilities can 
arise in a gravitating fluid because of thermal and lepton 
transport that would otherwise be stable in their absence. 
These instabilities, involving the transport of two quan- 
tities (e.g., energy and leptons), are referred to as doubly 
diffusive instabilities. 

To introduce the phenomenon, let us consider several 
classic examples of doubly diffusive instabilities. A dou- 
bly diffusive instability will arise if there are two differing 
"substances," say heat and leptons, with gradients that 
can act oppositely in their stabilizing effects. For the 
sake of a concrete model, consider the case of heat and 
salt in water. Let us first imagine the case in which a 
region of hot salty water overlies a region of cold fresh 
water (jSternl IT960I) . Normally the diffusion of heat oc- 
curs more rapidly than the diffusion of salt, and we will 
assume this to be the case here. Let us also assume 
that the magnitudes of the gradients in heat and salt 
are such that the fluid is stably stratified gravitationally 
(i.e., Ledoux stable). Imagine a parcel of hot salty water 
to be depressed slightly, as shown in Figure 1. The rapid 
diffusion of heat will thermally equilibrate this parcel of 
water with the background, but the slower diffusion of 
salt will result in the water remaining salty. We thus 
end up with a pocket of cold salty water in a background 
cold fresh water, and this pocket, being denser than the 
background, will continue to sink. The result is that "fin- 
gers" of salty water will poenetrate the fresh water on a 
thermal diffusion time scale. Note that what creates the 
instability in this otherwise stable situation is the pres- 
ence of diffusion, and thus its description as a "diffusive 
instability." The above described instability is referred 
to as salt fingers, and is an example of the kind of doubly 
diffusive instability that occurs if a slowly diffusing sub- 
stance with a destabilizing gradient is stabilized in the 
absence of diffusion by the gradient of a rapidly diffusing 
substance. 

Now imagine that the spatial configuration described 
above is reversed, namely, that a region of cold fresh wa- 
ter overlies a region of hot salty water, so that the gra- 
dient of salt is now stabilizing, the gradient of heat dsta- 
bilizing, and the magnitud es of the grad ients such that 
the fluid is Ledoux stable (|Veronislll965l) . Imagine fur- 
ther that a parcel of cold fresh water is depressed slightly 
into the hot salty water, as shown in Figure 2. Again the 
rapid diffusion of heat will thermally equilibrate the par- 
cel with the background, but the slower diffusion of salt 
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Fig. 1. — Simple model of the salt-finger instability. 



will result in the parcel remaining fresh. We thus end 
up with a parcel of hot fresh water in a background of 
hot salty water, and the parcel, being less dense than 
originally from the background, will accelerate upwards 
faster than would have been the case with no diffusion. 
This can result in an overstable situation, driven by the 
phase lag between the temperature and velocities of the 
fluid e lement which causes net work to be done over each 
cycle ( Spiege] Hl972ft . In this case the parcel will oscillate 
with growing amplitude. This instability is referred to as 
"semiconvection," and is an example of the kind of dou- 
bly diffusive instability that occurs if a rapidly diffusing 
substance with a destabilizing gradient is stabilized in 
the absence of diffusion by a slowly diffusing substance 
with a stabilizing gradient. 

Doubly diffusive instabilities such as described above 
can occur in the co llapsed core of a supernova progenitor 
()Smarr et alJl!98ll) . Neutron rich (low Y e , where Y e is 
the proton fraction) matter tends to be denser than neu- 
tron poor (high Y G ) matter at the same temperature and 
pressure because the latter contains more light particles 
(electrons and electron neutrinos) that contribute to the 
pressure. Thus neutron richness in stellar cores can be 
thought of as playing a role analogous to that of salt in 
the discussion above. The ramp-up of the bounce shock 
on core rebound results in an outwardly positive (stabiliz- 
ing) entropy gradient throughout much of the inner core. 
On the other hand, a negative lepton gradient formed at 
and initially confined to the vicinity of the zvsphere im- 
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Fig. 2. — Simple model of the semiconvective instability. 



mediately after shock breakout is advected deeper into 
the core with time because of the inward matter flow 
and outward diffusion of z/ e 's. A destabilizing negative 
Y e gradient stabilized by the outwardly positive entropy 
gradient can set the stage for a doubly diffusive instabil- 
ity analogous to salt fingers if lepton transport is slower 
than thermal transport. This instability is referred to as 
"neutron fingers." 

The Livermore group has reported successful explo- 
sions in their spherically symmetric supernova simu- 
lations employing sophisticated multi-energy neutrino 
transport, and have emphasized the importance of the 
neutron-finger insta bility in powering the s e exp losions 
HWilson et n ,UI198?i IWilson MavldllflSa ITflfll. The 
issue of whether or not this instability is operative, 
and if so whether or not it will power an explosion 
is important, as no other supernova simulations em- 
ploying neutrino transport of comparable sophistica- 
tion but witho ut invoking; neutron fingers has produced 
an explosion (|Bruennl f 19931: iM ezzacao pa et alJ 120011: 
iLiebendorfer et aJ.ll200H iBuras et alJl200fl . 

The Livermore group argue that thermal equilibration 
between a fluid element and the background is driven 
by the total neutrino number density (i.e., ^ e 's, z^'s, 
and zVs and their antiparticles), whereas equilibration 
in composition (i.e., lepton fraction, Yg) is driven only 
by the difference in the number densities of v^s and z/ e 's. 
The ratio, R, of the two equilibration rates is thus given 
approximately by 



(1) 



where rii is the number density of neutrinos of type i. 
From their simulations they typically find that R < 0.1. 
Based on the magnitude of this inequality the Livermore 
group tests for the neutron-finger instability by displac- 
ing a fluid element at constant composition but equi- 
librated with the background material in pressure and 
temperature. If the density of the fluid element rela- 
tive to its surroundings after displacement is such as to 
drive the displacement further, the fluid is taken to be 
neutron- finger unstable. Mathematically, their criterion 
for neutron-finger instability is given by 



dY e J rji p 



dr 



> 



(2) 



where is the Y e gradient of the background. If a 
region was established to be neutron-finger unstable, a 
mixing length algorithm was then employed to simulate 
the convective fluid motions and mixing. The Livermore 
group found that the result of the convective-like motions 
due to the neutron-finger instability in their supernova 
simulations was the production of a larger v e luminosity, 
as neutrinos and electrons previously trapped in the core 
were brought to the vicinity of the neutrinosphere. This 
increase in v e luminosity was found to be sufficient to 
power an explosion. 

In several p ap ers, iBruenn et alJ (|l995l) and 
iBruenn fc Dineval (|l996t ) questioned the assump- 
tion that the ratio, R, of the rates of composition to 
thermal equilibration is given by equation Q). Situa- 
tions were envisioned in which a fluid element perturbed 
with respect to the background could drive counter 
flows of ^ e 's and z/ e 's so that the numerator of equation 
JU would be additive in n Ue and n Pe and the denomi- 
nator subtractive . Deta iled simulations, performed by 
IBruenn fc Dineval (Jl99 6). of the equilibration of fluid 
elements perturbed with respect to the background 
showed that in many cases counter flows of z/ e 's and P e 's 
did indeed occur. Furthermore, it was found that lepton 
transport proceeded quickly by the rapid transport of 
low energy z/ e 's and zVs, as the opacity for these was 
relatively small. On the other hand, thermal equilibra- 
tion, being more efficiently accomplished by high energy 
neutrinos, was slowed by the relatively large opacities 
encountered by these neutrinos. Finally, the z^'s, ^ r 's, 
and their antiparticles were less effective at mediating 
thermal equilibration because of there weak thermal 
coupling with the matter. The result was that the ratio 
of composition to thermal equilibration depended on 
how the fluid element was perturbed with respect to the 
background, and in many cases this ratio was greater 
than unity. 

Given its potential role in the core collapse supernova 
mechanism, there is clearly a need to clarify the issue 
of the neutron-finger instability in a proto-supernova. 
Additionally there has been the advent during the 
past decade of mu lti-dimension a l cor e collapse su- 
perno v a modeling ([Miller et all Il993t iHerant et alJ 
[ 19941: iBurrows et al] 
iMezzacappa et alj|1998a 



vliller^5|Jl993; 
jay;|Fryerl|199fH 



Fryer WarrerJ 
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l2Q02t iBuras et alJ 12003). with more simulations of 
greater sophistication to come. The extant simulations 
show the presence of fluid instabilities both above 
and below the neutrinosphere. The fluid instabilities 
above the neutrinosphere are driven by neutrino heating 
immediately above the gain radius, which is the surface 
above the neutrinosphere where neutrino cooling and 
heating balance. This instability is analogous to that of 
a fluid in a pot heated from below. The fluid instabili- 
ties below the neutrinosphere are more complicated to 
characterize, as neutrino transport drives both thermal 
and composition changes and the fluid instabilities can 
b e doubly diffusive. 

iBruenn Dineval (|1996T) developed a procedure for 
characterizing the stability /instability of a fluid element 
as a function of the entropy and composition gradient 
of the background and found that for typical conditions 
inside the collapsed stellar core of a supernova progen- 
itor the fluid was either unstable to semiconvection or 
stable. However, they sampled only a small region of 
the run of thermodynamic conditions inside a supernova 
progenitor. In this paper we wil l furth er develop the 
methodology of IBruenn fc Dine"val (|l996ft and use it to 
examine the stability of proto-supernovae for the entire 
run of thermodynamic conditions from the center of the 
core to the neutrinosphere. Fut her more, for each set of 
thermodynamic conditions we will consider the stabil- 
ity as a function of the size of the fluid element. The 
latter is important in that it governs the rate of thermal 
and composition equilibration. In Section[2]we begin our 
discussion of fluid stability in the ambience of neutrino 
mediated thermal and lepton transport by presenting the 
equations of motion of a fluid element that we will use to 
analyze stability. These equations involve four response 
functions that describe how thermal and lepton diffusion 
is driven by an entropy or a lepton fraction difference 
between a fluid element and the background. In Sec- 
tion 03 we describe how the solutions of these equations 
of motion are used to define various types of instabil- 
ity. We also examine the stability of a fluid element as a 
function of the entropy and lepton gradients of the back- 
ground for various representative sets of the "response 
functions." These examples illustrate the different kinds 
of instabilities that can arise in a proto-supernova. In 
particular, we find that neutron fingers is unlikely to oc- 
cur anywhere below the neutrinosphere. However, we 
discover and describe several new and potentially impor- 
tant modes of doubly diffusive instabilities which we refer 
to as "lepto-entropy fingers" and "lepto-entropy semi- 
convection." Both of the instabilities involve the cross 
response functions in essential ways. In Section 0] we 
describe how the response functions are computed for a 
given thermodynamic state and fluid element size. In 
Section we begin the application of our stability anal- 
ysis to core collapse models by describing these models, 
and in Section|£|we present the results of an extensive set 
of calculations giving the stability or type of instability 
and its growth rate as a function both the fluid element 
size and its location in the core. In Section [3 we com- 
pare our analyses and results with prior work, and in in 
Section [HI we present our conclusions. 



We consider a mix of hot, dense material and neutri- 
nos at conditions typical of the region below the neutri- 
nosphere of a proro-supernova (p > 10 12 g cm -3 , T > 10 
MeV, Yi < 0.5. Such a mix can be considered to be 
equilibrated with respect to strong, electromagnetic, and 
weak interactions. Its thermodynamic and compositional 
states can therefore be specified by three independent 
variables, which we take to be the pressure p, entropy s, 
and lepton fraction fraction Y. 

In a manner analogous to iGrossman et all (|l993t ). we 
consider the behavior of an individual fluid element in 
a background which is assumed to be time-independent 
and one-dimensional. The z-axis is taken normal to the 
plane of symmetry and the background is assumed to 
have constant gradients in entropy s and lepton fraction 
Yg. We will assume that the fluid element is always in 
pressure equilibrium with the background, i.e., p = p, 
where unmarked and barred variables refer to the fluid 
element and the background, respectively. However, the 
fluid element and the background may differ locally in 
s and Yg, and these differences will be denoted by the 
quantities S and 0y g , where S = s (z) = s(z) — s(z) 
and 0y e = 6y e (z) = Yg(z) — Yi(z). Given p(z), s(z), and 
Yi(z\ the thermodynamic state of the fluid element is 
given by its values of S , 0y g1 and its vertical location z 
(positive for directions away from the center of the core) 
through 



s(zq) + 



s(z) = s(z) 
ds 



dz 



■(*) 

(z-z ) + e s (z) 



(3) 



Y i (z) = Y e (z) + Ye (z) 

(z-z )+9 Ye (z) (4) 



p = p(z). 



(5) 



The quantities S1 and 0y £ , will change with time by 
the effect of energy and lepton transport by neutrinos, 
and by the vertical motion of the fluid element through 
the gradients of s and Yg of the background. We assume 
that neutrino energy and lepton transport is linear for 
small # s , and 0y £ and describe it by the four "response 
functions" E s , Ey £ , T s , and Ty £ , defined by 



0Yo 



0y i 
0y £ 



e s =o 



o s =o 



(6) 



Thus, E s gives the rate of change of S due to the pres- 
ence of a nonzero entropy difference, S , and a zero lepton 
fraction difference, 0y £ , between the fluid element and the 
background, T l y i gives the rate of change of 6 S due to the 
presence of a nonzero lepton fraction difference, 0y n and 
a zero entropy difference, S , and so forth. With these 
definitions, it follows that 
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and 



where 



Y,VY, 



(8) 



and 6y e are the rate of change respectively 

V 

of S and 0y e due to neutrino transport. The presence 
of the "cross" response functions Ey^ and T s is due to 
the fact that, (1) entropy transport is brought about by 
both thermal and lepton transport, (2) that neutrinos 
transport both energy and leptons, and (3) that trans- 
port is inherently a dissipative process. Thus, for exam- 
ple, a difference, 6y e , in Yg can force a rate of change 
of s as well as Yg, as can a difference, (9 S , in s. These 
cross response functions introduce important additional 
patterns of stability /instability behavior as a function of 
the gradients in s and Y^, as will be described below. Dis- 
cussions of doubly diffusive instabiliti es apart from the 
proto-supernova problem and prior to iBruenn fc Dineval 
(1996) did not relate to a context in which cross response 
functions were important, and these were not considered. 

The equations for S and 0y e are completed by adding 
to equations J7J) and J5J) the effect of the vertical motion 
of the fluid element, v = i, to get 



6 S — Hi, s 



J Ye t 



ds 
dz 



r Ye 0y e -—v 



(9) 



(10) 



The equation of motion for the fluid element is given 
by the buoyancy acceleration driven by the difference 
between its density and that of the background, viz., 



dp 

~d~ s 



dp_\ 



(11) 



where g is the magnitude of the acceleration of gravity. 

3. STABILITY ANALYSIS 

To investigate the stability of a fluid element described 
by equations (0 - , we regard these equations as giv- 
ing the time evolution of 6 S , 6y £ , and v and consider in 
this Section classes of possible time dependent behav- 
iors of a fluid element perturbed from equilibrium as a 
function of the gradients in s and Yg, the radius, i?biob> 
of the fluid element, and the gravitational force, g. In 
particular, we look for solutions of the form 



= sO e 



at 



= YiO e a , v = v e 



at 



(12) 



where the eigenvalues, a, will determine the stability of 
the fluid, as descr i bed below. 

In Sections 13.11 - 13.31 we look at some particular cases 
where the response functions take on simplified and spe- 
cial sets of values. The purpose is to relate some clas- 
sic examples of instability to the solutions of equations 
equations (0 - ([TT|). and to build from the simple to the 
complex. We delay until Section 13.41 an analysis of the 
general case in which all four response functions have ar- 
bitrary nonzero values. The reader who is interested in 
proceeding directl y the the general case is advised to skip 
directly to Section 1^1 



3.1. Ledoux Stability Limit 

We begin with a well-known case by considering the 
limit in which there is no thermal or lepton transport 
(i.e., the limit in which all four response functions are 
zero). Instability in this case is referred to as Ledoux 
instability, and equations ® - (fTT]) in this case reduce to 



ds 
dz 

dz 



dp_ 

dY e 



lY e - 



(13) 
(14) 
(15) 



The eigenvalues a are given by 



a 






a 



(dp\ g(dp\ 
\ ds Jp,Yt p V ay W 



(16) 



with the solutions 



01.2 ±V^4 



'din p 



d\nl 



dhis) vYi dz 



J Y,i 



din p 

dhiYi 

<J3 = 0. 



dlnYi 

dz 

(17) 



where we define uo s and uoy i by 



-9 



'<91np^ dins 

, din s . 



dz 



din p \ dlnYi 



K d\nY £ 



p,s 



dz 



(18) 



Physically, if ^p- = then uj s is the frequency that a 
fluid element would oscillate about a stable equilibrium if 
uol > 0, and the growth rate that would characterize the 
motion of a fluid element away from an unstable equilib- 
rium if uol < 0. coy £ has an analogous interpretation, e.g., 
if = then uiy e is the frequency that a fluid element 
would oscillate about a stable equilibrium if co Yi > 0? 
and the growth rate that would characterize the motion 
of a fluid element away from an unstable equilibrium if 
u Y < 0. u s and uoy g thus characterize the effect of the 

background gradients, and d , respectively, on 

the dynamics of a fluid element perturbed from equi- 
librium. In equation ([T7)l. the solution as = arises 
from the fact that a fluid element at rest and having the 
same thermodynamic conditions as the background (i.e., 
6s = 6y £ = v = 0) will remain at rest. Note, however, 
that the equilibrium implied by the solution as = may 
not be a stable equilibrium. 

The criterion for Ledoux stability is that the eigenval- 
ues a be imaginary, which requires the following inequal- 
ity be satisfied: 



din p 
dins 



din i 

dz 



dlnp 
dlnYo 



dlnYi 

dz 



< 0. (19) 
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If this inequality is satisfied, then the two solutions, 6J1.2, 
can be written 



(71.2 = ±id)B-V 



(20) 



where 



ub-v 




is the Brunt- Vaisala frequency. In this case the fluid el- 
ement will oscillate about the equilibrium configuration 
with the angular frequency ojb-v and constant ampli- 
tude. This situation is therefore stable. On the other 
hand, if inequality (fTH|) is not satisfied, i.e., if 



<91n p 
din s 



dlnl 

dz 



din p 
din Y £ 



dlnYf 



> 0, (22) 

p,y £ az \ ULLL1 t/ P ,s 

then the solutions, G\ and <J2, can be written g\^i = 
±\oub-v\-> with the positive solution corresponding to 
an exponentially growing amplitude with growth rate 
\u B -v\ given by 



\ub-v\ 



= V9\ 



'din p 
t dins 



dins 

dz 



din p 
dlnYg 



p,s 



dlnYg 

dz 



.(23) 



In this case the fluid is said to be convectively unstable. 
Thus, in the case of no thermal or lepton transport, the 
"Ledoux critical line" 



din p 
dins 



dins 



din p 
dlnYg 



dlnYg 

dz 



(24) 



separates, in the - ^ - plane, a convectively stable 
region where the left-hand side of equation (|24[) is < 
from a convectively unstable region where the left-hand 
side of equation (|2^|) is > 0. 

It is interesting to consider the signs of the coeffi- 
cients of ^ and in equation The derivative 
/ ainp \ • equal to — Ts— , where c v = T f^fr) is 

the always positive constant pressure specific heat and 
a = \(^r) pYg * s the expansivity. Provided the lat- 
ter is positive, which is almost always the case, the 
derivative ( ) is negative, and a positive gradi- 
ent is s is therefore stabilizing. On the other hand, the 
derivative ( q^y ) is equal to —Ygn s ( Jy- ) , where 

is the isentropic compressibility (al- 



v y op 



ways positive), so the sign of ( d\nY e ) can ^ e ne g a ^ive 

or positive depending on the sign of (jjy^ • The latter 

can be positive or negative depending on the magnitude 
of Yg. To explain this curious behavior, note that in- 
creasing Yg increases the number of light particles while 



leaving the baryon number essentially unchanged. The 
entropy, which is being held constant in (jyj) > de- 
pends on both the number of particles and the temper- 
ature. Increasing the former will typically (but not al- 
ways) require a decrease in the latter in order to maintain 
constant entropy, and the pressure can therefore increase 
or decrease. Whether p increases or decreases depends on 
which particle specie is dominating the pressure. Figure 
03 shows the contribution of the various particle species 

to . It is seen that the light particles (e 's and 

V 1 J p,s 

zVs) provide a positive contribution to ( Jp- ) for all 

values of Yg for the given thermodynamic conditions, due 
to their partial degeneracy, but the baryons are nonde- 
generate and provide a negative contribution. For values 
of Yg above 0.11, and for the indicated values of p and s, 

the contribution of e _, s and zVs to ( 4& J dominate 



\dY t J 



and this derivative is positive. But for values of Yg below 
0.11 the baryons dominate and this derivative is nega- 
tive. For values of Yg close to 0.11 this derivative, and 

therefore ( £ lrL y ) , is almost zero. The stability of a 

V n 1 J p,s 

fluid element in this latter regime is extremely insensi- 
tive to the gradient dl ^ Yl , and the Ledoux critical line is 
almost vertical. 



0.8 
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Fig. 3. — Contributions to the derivative 



The dashed 



lines show the contributions of the indicated particles to this deriva- 
tive, the solid red line shows this derivative with all contributions, 
and the solid blue line shows the corresponding values of Y e . 



To exemplify the above discussion, Figure SI shows the 
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Fig. 5. — Ledoux stability lines for s = 4 and the indicated 
values of Yg. 



consider the case in which the former are nonzero while 
the latter are still zero. (That is, we take E s 7^ 0, Ty e 7^ 
0, Yiy £ = 0, and T s = 0.) This corresponds to the case 
in which an entropy gradient drives an entropy, but not 
a lepton flow, and a lepton gradient drives a lepton but 
not an entropy flow. Then equations Q - JUJ become 



_ . ds 
Z S 6 S - —v 
dz 



^ = _9_ (dp 

p yds 



T Yi Yi - —v 



9_(dp_\ 



P,Y £ 



P \9Y £ J 



6y 



(25) 



(26) 



(27) 



Let us further suppose in this Section that Ty e — > — 00, 
that is, that the fluid element is instantaneously equili- 
brated in Yg with the background, leaving for the next 
Section the case in which Ty £ is finite. If Ty £ — >• —00, 
it follows that 0y e — >• which implies that dy t — > 0, so 
equation must reduce to 



= T Yi 0y a - — 1 



(28) 



dYi 



i.e., 6y g — > and Ty £ — >• — 00 such that Ty^y £ — >• ^g^. 
The equations for the nonzero variables now simplify to 



(22; 



ds h,Ye 



(29) 



(30) 



Using the first and the last of equations (|T2|) . the condi- 
tion for nontrivial solutions is 







<91n p 
<91n s 



dz 



(31) 



Let us consider the case E s 
If E s = 0, then 



and £ s 7^ separately. 



Ledoux critical line as a function of Yg for the indicated 
values of p and s. As Yg decreases it is seen that the slope 
of the Ledoux critical line goes from negative through the 
vertical to positive. The meaning of a positive slope for 
the critical line is that a negative (destabilizing) gradient 
in s is now balanced by a negative rather than a positive 
gradient in Yg. At lower entropies, on the other hand, 
such as the case of s — 4 shown in Figure the baryons 

do not dominate the derivative (jyj) even for Yg as 

low as 0.07, and the slope of the Ledoux critical line 
remains negative for all values of Yg shown. 

The sign change in the slope of the Ledoux critical 
line as Yg falls below a critical value will have important 
ramifications below for the locations of the instability 
regions in the - d plane when we consider more 
general sets of values for the response functions. 

3.2. Schwarzschild Stability Limit 

Referring to E s and Ty £ and to Ey € and T s as the 
direct and cross response functions, respectively, let us 




± 



\ 



-9 



dlnp 
din s 



d\ns 

dz 



(32) 



where in the last expression we have explicitly taken into 
account the fact that (fr^) < 0. Equation 

expresses the Schwarzschild stability criterion, viz., if 
> then the roots are imaginary and the fluid is 
oscillatory stable, the oscillation frequency being u s de- 
fined in equation (|T8|) ; if ^| < then there is one positive 
and one negative real root and the fluid is convectively 
unstable with a growth rate, \lj s \, given by the magnitude 
of the positive root, e.g., 



\ 



<91n p 
d\n s 



dins 



dz 



(33) 



8 



The critical stability line, i.e., the line in the ^f^-^ 1 
plane that separates real from imaginary roots (instabil- 
ity from stability) is the vertical line = 0. 

Assume now that Ty £ is as above (viz., Ty £ — > — oo) 
but E s 7^ 0, and that E s < 0, i.e., that neutrino transport 
will drive an entropy perturbation of a fluid element back 
towards equilibrium with the background. In this case 
the solutions of equations (|25l - (|27[) are 



Close but to the left of critical line 1 separat- 
ing stable from unstable solutions, where |^p| << 



(dlnp\ 
\ dlnS Jp,Yi 



, or equivalently where \uj s \ « E s /2, 



the growth rate, \ws\, of the instability is 



\vs\ 



(36) 



E* 
2 



db 



51 
4 



db 



51 
4 



dlnp 
dins 



P,Yt 



dins 

dz 



(34) 
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Fig. 6. — The two critical lines characterizing Scwarzschild sta- 
bility. 



In contrast to the case in which E s = and for which the 
critical stability line = separates both the real and 
imaginary solutions and the unstable and stable regimes, 
with E s < the critical line bifurcates into two paral- 
lel critical lines. The situation is shown schematically in 
Figure El The critical line = 0, which we will refer 
to as critical line 1, is the line along which one real root 
is zero, and so changes sign as this line is crossed. It still 

separates stable from unstable solutions. However a sec- 

-l 

ond critical line = -j^ 

dz 4g 



now appears, 



which we refer to as critical line 2. Along this line two 
real roots switch to a complex conjugate pair, so this crit- 
ical line separates real from complex solutions. To the 
right of the critical line 2 the fluid is stable; if perturbed it 
will oscillate with decreasing amplitude. Along critical 2 
line the neutrino transport critically damps the fluid. In 
the region between the two critical lines the fluid is sta- 
ble but non-oscillatory — neutrino transport overdamps 
the fluid. To the left (^p < 0) of the critical line 1 the 



fluid is unstable, as in the E, 
can be written 



= case, and equation 



2 V 4 



(35) 



The growth rate is thus reduced from the E s = case 
by a factor of . This factor arises from the reduction 
in S (and thus in the buoyancy force) due to the ten- 
dency of neutrino transport (as modeled by a negative 
E s ) to thermally equilibrate the fluid element with the 
background. Finally, far to the left of the first critical 



I d In a 
I dz 



r 2 

» A~ 



fdhip\ 



line, where 

where \u 3 \ » E s /2, the growth rate of the instability is 



Yi 



, or equivalently 



us 



Es 
2 ' 



(37) 



i.e., just slightly less than the E s = case. Here there is 
also a reduction in 6 S (and thus in the buoyancy force) 
due to the tendency of neutrino transport to thermally 
equilibrate the fluid element with the background, but 
this effect is now small compared with the growth rate, 
\w 8 \, the fluid element would experience in the absence 
of transport. 

Now imagine for a the moment that E s > 0, i.e., that 
neutrino transport will drive an entropy perturbation of 
a fluid element further away from equilibrium. Then 



2 



db 



\ 



din p 
din s 



P,Yt 



dins 

dz 



2 



db 



51 
4 



(38) 



In this case there is no stability. Critical line 1 (^^- = 0) 
separates the region to its left where equation (pT5j) gives 
one positive and one negative root from the region to 
the right but to the left of critical line 2 (^1 = 

-l 

) where equation gives two posi- 



( dlnp\ 
\ dlns Jp,Yi 



4g 

tive roots. To the right of critical line 2 equation ((38)) 
gives complex solutions with a positive real part — the 
motion of a perturbed fluid element in this case is oscil- 
latory with a growing amplitude, i.e., semiconvective. 

3.3. E s ^0,T Y£ ^0 

We generalize the case considered in the preceding Sec- 
tion by allowing Ty £ to be finite rather than - oo; that is 
we consider the case in which the two "direct" response 
functions, E s and Ty^, are nonzero and arbitrary, but 
the two "cross" response functions, Ey^ and T s , are still 
zero. Then equations (0 - become equations (|25[) - 
(|27)l . which are written here again as 



6 S — £c0c 



ds 

dz 



(39) 



9 



and 



-{¥) 9 >- g -(w) 



(40) 



(41) 



Solutions have the form given by equations ifT2jl . and 
exist for values of a satisfying 



where 



a 3 + Aa 2 + Ba 1 + C = 0, 
A=-(X s + T Ye ), 



(42) 
(43) 



l(AB -3C)- l~A 2 . 
6 V ; 27 



(48) 



Thus, critical line 2 is given by equations - 

(|47j) and (|48j) . This is a fairly messy nonlinear equation 

in Yl and and will not be written down here. 

dz dz 

It will be solved numerically for a number of illustrative 
examples below. 

Finally, and unlike the case of Schwarzschild stability, 
there is the possibility that the real part of a complex 
conjugate pair of roots can change sign. This will hap- 
pen if there is a complex conjugate pair of roots, and if 
they sum to zero. To determine the condition for this 
Stegunlll965D that the roots of cu- 



note (|Abramowitz 
bic equation (|42|) satisfy 



B = Z s T Ye -g\ 
d\n p 



'dlnp\ dins 



-g 



\ dlnY £ 



dz 



(44) 



and 



C = Z s g 



dlnYf 



din p 
dlnY £ 

dlnp\ dins 



dins 



(45) 



The eigenvalues, cr, in this case are given by the solu- 
tion of a cubic equation, and so display a richer variety. 
In particular, there are now three roots rather than two. 
The three roots are either all real, or consist of one real 
root and a complex conjugate pair. In the above dis- 
cussion of Schwarzschild stability, critical line 1 in the 
d — plane was defined as the line along which 
a root is zero, so that a root changes sign as this line is 
crossed (i.e., a stable/unstable mode switches to a un- 
stable/stable mode). The analog of that line here is ob- 
tained by setting C in equation (|45|) to zero. This gives 



dlnY £ 

dz 



( d\np \ 
{dins) 



dins 



V din Y £ J 



(46) 



which we will also refer to as critical line 1. As in the 
above case of Schwarzschild stability, critical line 1 here 
is the line across which a root changes sign. 
The second critical line in the dl ^ Yf - - plane in- 

dz dz ^ 

troduced in our discussion of Schwarzschild stability was 
the line across which a pair of real roots switch to a com- 
plex conjugate pair (growing/decaying modes switch to 
growing/decaying oscillatory modes). The analog here 
is the line across which two real roots switch to a com- 
plex conjugate pair, and is o btained from the equation 
(lAbrarnowitz fc StemmllT965h 



q 3 + r 2 = 0, 



(47) 



(49) 



cri + cr 2 + cr 3 = -A, 
(J\(J2 + C r lC r 3 + ^"2^3 = B, 

Then = AB - C gives 

= aia 2 (cri + cr 2 ) + cricr 3 (cri + a 3 ) 

+^2^3(^2 + &3) + 2(71(72(73, 
= (<7i + Cr 2 )(<7i«7 2 + «7i«7 3 + <7 2 Cr 3 + (if ) 

where the subscripts 1, 2, and 3 can be cyclically per- 
muted in the above expressions. Thus equation (|KU|) is 
satisfied if any two roots sum to zero. This is the de- 
sired condition. Using equations - (|I5|) equation 
= AB — C can be written 



(50) 



( d\np \ 
\d\ns J 



P,Y e 



dins E a (E a + T y J 



d\nY £ 

g 

and we will refer to this line as critical line 3. 



( dlnp \ dz n ( d\np \ 



(51) 



10 



' c" I N p 





; 1 1 1 1 1 1 1 1 1 1 

V 


i i i | 
P = 


i , 1 1 1 1 1 1 1 , 1 1 1 

3x10 12 g cm" 3 _ 




\ \ NF \ 


s = 


5 






y\ = 


0.3925 


; 3 














Stable 




Convective 






; 


-10 2 s- 1 






- E yi = 







Semi-C 2^ 


- Y s = 









; Y yi = 


-3x1 2 s" 1 




\\l_edoux 


: 10 3 


3- 1 <|co BV |<10 4 S" 1 

, , , i , , , i , , , i , , 


, , 1 ,\ 


1 1 I 1 1 1 IX 1 *i« I 1 1 1 " 



2 

dins 

le "dz" 



Fig. 7. — Fluid stability as a function of the Ins and lnY^ 
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a few hundred milliseconds after core bounce. 
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thermodynamic state and values of S s , Ty £ , and ujbv- The value 
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a few hundred milliseconds after core bounce. 



the NF case, but the region resides on the unstable side 
of the Ledoux critical line and would therefore be unsta- 
ble even in the absence of thermal and lepton diffusion. 
We summarize our stability /instability classification in 
Table □ 

To characterize in more detail the stability /instability 
of the fluid in the dh ^ 1 - ^ff- plane for the conditions 
shown in FigurefTJ and to delineate the role of the critical 
stability lines, we will start in the upper right quadrant 
of that figure where the gradients of In s and In Yi are 
both positive. The condition that the fluid is stable, i.e., 
that all the roots of equation (|42|) are negative, if real, 
or that the real root is negative and the real part of a 
complex conjugate pair is also nega tive, is that all of 
the following inequalities be fulfilled ([Aleksandrov et all 
H96l : A > 0, C > 0, and AB > C, where A, B, and C 
are given by equations (f!5|) - (J35J). A, given by equation 
is clearly positive, as both S s and Ty e are negative. 
For the given thermodynamic conditions the logarithmic 

derivatives ( ) and ( J are both negative, 

therefore C, given by equation ((15)). is also positive. Fi- 
nally, 



Figures Q and |H1 show these three critical lines along 
with the Ledoux line in the - plane for the 

dz dz ^ 

indicated thermodynamic state and values of E s , Ty £ , 
and IcjbvI- The absolute value of the Brunt- Vaisala fre- 
quency, IcjbvI, where ojbv is given by equation (jSTJl. is 
an indicator of the buoyancy growth rate in the absence 
of neutrino transport. Its value, however, depends on the 
gradients of In s and In Yi and therefore varies over the 
d — ^jf- plane. A range of typical values for ubv is 
given in the figures. A comparison of these values with 
the absolute values, |E S | and |TyJ, of the direct response 
functions will indicate the relative importance of buoy- 
ancy versus neutrino transport on the dynamics of the 
fluid element. 



AB-C 



We have divided the 



din Y, 



plane into four dis- 
tinct regions delimiting different types of fluid stabil- 
ity/instability, the boundary between any two adjacent 
regions being one of the critical lines. In the region de- 
noted "stable" all roots of equation are negative, 
if real, or have a negative real part if a complex conju- 
gate pair. A perturbed fluid element in this region will 
return to its unperturbed state, and the fluid is there- 
fore stable. In the region marked "Semi-C" the largest 
growth rate arises from the positive real part of a com- 
plex conjugate pair of roots, and the dominate unstable 
mode is therefore oscillatory and growing. In this region 
the fluid is unstable to semiconvection. In the region 
marked "NF" the largest growth rate arises from a pos- 
itive real root, but the region resides on the stable side 
of the Ledoux critical line and would therefore be stable 
in the absence of thermal and lepton transport. The in- 
stability in this region is therefore driven by thermal and 
lepton transport. We refer to this instability for the mo- 
ment as "neutron fingers," although we will show below 
that for the conditions in the collapsed stellar cores of 
supernova progenitors that we have analyzed this insta- 
bility is not neutron fingers per se but something quite 
different. Finally, in the region marked "convection" the 
largest growth rate arises from a positive real root, as in 



^s9 



dhip 
dhis 



P,Y t 



d\ns 

dz 



din p 

iMYt 



dhiYt 
dz 



> 



(52) 



since each term to the left of the inequality sign of equa- 
tion ((52)) is positive. Thus all three of the above inequal- 
ities are satisfied in the upper right quadrant of Figure 
Q where the gradients of In s and In Y$ are both posi- 
tive, and the fluid is therefore stable in this region. For 
the particular case shown in Figure in the upper right- 
hand quadrant there is a negative real root and a complex 
conjugate pair with a negative real part. 

To characterize the stability /instability of the fluid in 
the remaining parts of the dl ^ e - ^jf- plane, let us 
move clockwise around the figure from the upper right- 
hand quadrant. The character of the roots (and therefore 
the stability of the fluid) does not change until critical 
line 3 is crossed, at which point the real part of the com- 
plex conjugate pair of roots changes sign from negative 
to positive. The fluid, if perturbed, will now oscillate 
with growing amplitude, and according to our taxon- 
omy is therefore unstable to semiconvection. Continuing 
clockwise, critical line 2 is crossed and the complex con- 
jugate pair of roots with positive real part become two 
real positive roots. The fluid, if perturbed, will move 
away from equilibrium with exponentially increasing am- 
plitude. The fluid is therefore unstable in this region, 
and since it is would be Ledoux unstable in the absence 
of transport, it is therefore convective. Crossing line 1 
at the bottom of the figure, one of the positive real roots 
changes sign, so there are now one positive real root and 
two negative real roots. Because of the positive real root 
the fluid is again convective. Continuing clockwise to the 
left of the figure and up, critical line 3 is crossed and the 
real part of a complex conjugate pair of roots, if present, 
would change sign. There not being a complex conju- 
gate pair of roots in this region, there is no change in 
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TABLE 1 
Stability Classification 



Classification 


Roots 


Subsidiary Conditions 


Stable 
Semi-C 
NF 
Convection 


ri < 
n > 0, si / 
n > 0, si = 

ri > 


Ledoux stable in the absence of transport 
Ledoux unstable in the absence of transport 



7*1 is the largest real root, or the real part of a complex conjugate pair, si is the imaginary part of the root of which r\ is the real part. 



the character of the roots and the fluid remains convec- 
tive. Continuing clockwise and crossing critical line 2 at 
the upper left-hand side of the figure, two negative real 
roots change to a complex conjugate pair with a negative 
real part, leaving one real positive root. The fluid is still 
convective. On crossing the Ledoux line the character of 
the roots remains unchanged but the fluid would be sta- 
ble in the absence of thermal and lepton transport. The 
fluid is now unstable to neutron fingers. Finally, crossing 
critical line 1 from left to right at the top of the figure 
the positive real root changes sign and we are left with a 
negative real root and a complex conjugate pair with a 
negative real part. We are back in the region where the 
fluid is stable. 
The above taxonomical sectoring of the dl ^ Yf - 

_ & dz dz 

plane in Figure [7| into regions of stability, convection, 
semiconvection, and neutrino fingers is a general charac- 
terization of a gravitating fluid with thermal and lepton 
transport represented by negative values for E s and Ty e . 
The location of these regions will depend, however, on 
the magnitudes of S s and Ty £ , and on the thermody- 
namic state of the fluid and the force of gravity. Figure 
[7| has illustrated a case in which Ty e = 3H S (lepton equi- 
libration more rapid than thermal equilibration). Inter- 
changing the values of E s and Ty £ , so that Ty e = 
(thermal equilibration more rapid than lepton equilibra- 
tion), and keeping the same thermodynamic state and 
gravity gives regions of stability and instability shown 
in Figure |H1 Note, as is apparent from equations ((46)) 
and (pTTJ). that interchanging the values of E s and Ty £ 
interchanges the slopes of critical lines 1 and 3 (and also 
changes the intercept of critical line 3). Also, the lo- 
cation of the NF and Semi-C regions are interchanged. 
Referring back to the discussion in Section [I] Figures [7| 
and [SI show, respectively, that a destabilizing (negative) 
gradient in Y$ stabilized by a stable (positive) gradient in 
s leads to semiconvection if Ty £ > S s (lepton transport 
more rapid than thermal transport) and neutron fingers 
if E s > Ty £ (thermal transport more rapid than lepton 
transport). The opposite is true if the destabilizing (neg- 
ative) gradient is s and the stabilizing (positive) gradient 
is Yi (tops of Figures Q and (HJ). 

Figures 03 and E| show the effect of successively increas- 
ing the values of E s and Ty £ , but keeping their ratio the 
same, and keeping the thermodynamic state and gravita- 
tional acceleration the same (all equal to the case shown 
in Figure [7J). It is apparent from equations ((46)) and (|51)) 
that the slopes of critical lines 1 and 3 are unaffected 
by the magnitudes of S s and Ty £ provided that their 
ratio remains the same. However, the intercept of crit- 
ical line 3 with the vertical axis becomes more negative 
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dients for the case in which |TyJ = 3|S S |, but larger in magnitude 
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(recall that ^ d\nY £ ) ^ s ne g a ^ive for the case being con- 
sidered). The latter effect is to extend the stable region 
with increasing magnitudes of E s and Ty e into the region 
that was formerly semiconvective. 

An effect not evident from the figures is the dependence 
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of the instability growth rates on the magnitudes of S s 
and Ty £ for given values of the gradients dlnY £ /dz and 
d In s/dz. Qualitativelly, the growth rates for instabilities 
that are driven by thermal or lepton transport (semicon- 
vection or neutron fingers) increase with the magnitudes 
of E s and Yy £ , while those that are buoyancy driven 
(convection) decrease with the magnitudes of E s and 
Ty £ . The latter effect arises because thermal and lepton 
transport tend to equilibrate the perturbed fluid element 
with its surroundings, thus reducing the buoyancy forces 
driving the instability. More quantitative results for in- 
stability growth rates will be presented for the surveys 
of proto-supernovae in Sectional 
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Fig. 11. — Fluid stability as a function of the Ins and In gra- 
dients for the case in which E s and Ty € have the same magnitudes 
as in the case shown in FigureQbut opposite (positive) signs. The 
thermodynamic state and values of ujbv are the same as in the 
case shown in Figure 171 

Finally, we show in Figure JTTJ the effect on fluid sta- 
bility of positive values of S s and Yy £ . In this figure, 
the thermodynamic state and gravitational acceleration 
are the same as that shown in Figure J7J), but both the 
values of E s and Ty £ have the same magnitudes but pos- 
itive signs. This causes the fluid to be destabilized for 
all values of dlnYg/dz and din s/dz, and is an alogous to 
the case discussed at the end of Section 13.21 While we 
have not encountered positive values of S s and Ty e in 
the proto-supernovae we have surveyed, we do encounter 
positive values of the cross response functions T i y i and 
T s , as will be discussed in the next Section. 

3.4. General Case - Lepto-Entropy Fingers and 
Lepto- Entropy Semiconvection 

In the general case all of the response functions, S s , 
Ey £ , T s , and Yy £ , given by equations (0) are nonzero, 
and the equations describing the motion of a perturbed 
fluid element are equations (0 - GU- Solutions for the 
motion of a perturbed fluid element again have the form 
given by equations (|T2|) . and exist for values of a satis- 
fying 



a 3 + Act 2 + Ba 1 + C = 0, 
where A, £?, and C are now given by 



(53) 



(E. + Ty £ ), 



B 



din p 
din s 



dins 



Sy £ T s 
din p 
d\nYo 



d\nY e 



(54) 



, (55) 



~ Ts9 Y f 



din p 
dlnY £ 

din p 
din s 

din p 
din s 

din p 
dlnYp 



p,s 



dlnY £ 

dz 



dim 



P,Yt 



p,s 



dz 

dlnY £ 

dz 

dins 

dz 



(56) 



Critical lines 1 (across which a root changes sign) and 3 
(across which the real part of a complex conjugate pair 
changes sign) are now given, respectively, by 



dlnYj 

dz 



<91n p 
din Yg 



( d^p\ 
\dlns ) 



P,Yi 



dins 



^Y, Y£ ( §u) - E s ( ) 



dz 



(57) 



and 



dlnY £ 

dz 



dlnp 
din s 



P,Y £ 



s y ( dlnp \ 
Y t L « \d\nY e ) ps dim 



(alna) p . 



1 (T s £y f -T F ,S S )(S S 



\ din Y e J 



Q Yi^ ( dlnp \ 1 y ( gin p \ 



(58) 



Nonzero values of the cross response functions, T i y i and 
T s , if no n- negligible, can significantly affect the existence 
and location of the stable region, and the location and 
nature of the three unstable regions in the ^r-^- - dln Yl 



dz 



plane. Figures El - El show cases in which the thermody- 
namic state, the value of the gravitational acceleration, 
and the two direct response functions, S s and Ty £ , are 
the same as for the case shown in Figure but the cross 
response functions have various nonzero values. The rel- 
atively small negative values given to Ey^ and T s in Fig- 
ure El cause the regions of stability and instability to be 
substantially shifted relative to that shown in Figure 
Changing the sign of the cross response functions from 
negative to positive causes further shifts in the regions 
of stability and instability, as evident from a comparison 
of Figures d and C21 

Increasing the value of T s by an order of magnitude, 
which is more consistent with the numerical results for 
the response functions, leads to some interesting results. 
The stable region collapses to a very small sector about 
vertical din s/dz axis, as shown by Figure lT^l and several 
new modes of instability appear. Figure IT5l shows a case 
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Fig. 12. — Fluid stability as a function of the Ins and lnYg 
gradients for the case shown in Figure^] but with Sy £ = —5 x 10 2 
s- 1 and T s = -5 x 10 2 s -1 . 



Fig. 14. — Fluid stability as a function of the Ins and lnYg 
gradients for the case shown in Figure |§1 but with Ey £ = 5 x 10 2 
s- 1 and T s = 7 x 10 3 s -1 . 
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Fig. 13. — Fluid stability as a function of the Ins and lnY^ 
gradients for the case shown in Figure but with Sy £ = 5 x 10 2 
s- 1 and T s = 5 x 10 2 s -1 . 



for which the values of the response functions and the 
gravitational acceleration are the same as shown in Fig- 
ure EI but the thermodynamic state causes the deriva- 

tive ( jfi^y^ J to be positive rather than negative. (This 

can arise for small values of Yg, as discussed in Section 
13.11 ) It is seen that in this case the region of stability 
is reflected about the origin, and that again new modes 
of instability appear. As will be described below, these 
instabilities arise because of the substantial magnitudes 
of the cross response functions together the large ratio 
of T s to Sy^. These cross response functions cause a 
6y e (a difference in Yg between the fluid element and the 
background) to develop primarily in response to a S (a 
difference in s between the fluid element and the back- 
ground). This is in contradistinction to the Ledoux case 
in which a 0y £ develops from a displacement of the fluid 
element through a gradient in the background Yg. 

3.4.1. Lepto-Entropy Fingers (L-E-F) 



Fig. 15. — Fluid stability as a function of the Ins and lnYg 
gradients for the values of the grav itat ional acceleration and the 
response functions shown in Figure ITTI but for a thermodynamic 

state such that ( ^Sr) > 0. 



, din Y t , 



We will discuss the instabilities that that are exhibited 
in Figure fTTTI here and in the next Section, and return to 
the instabilities exhibited in Figure ITU in Section 13. 4. 31 
The instability on the right-hand side of Figure El de- 
noted by L-E-F is new in that it has not been described 
or analyzed in the literature up to now. Moreover, it is 
likely to play the dominate role in producing convective- 
like fluid motions below the neutrinosphere of a proto- 
supernova. We begin by noting that the relative values 
of the response functions chosen for Figure IT51 as well as 

the fact that ( o\nY e ) > ^ are re P resen tative of these 

parameters in the region below the neutrinosphere and 
above the cold, unshocked inner core of a supernova pro- 
genitor. The relative values of the response functions 
seem to be generic, and will be discussed in more de- 
tail below. The positive value of ( ) is due to 
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low values of Y$ which come about because the mate- 
rial in this region, unlike that in the cold inner core, 
encountered the shock and did so while still outside the 
neutrinosphere. This material therefore suffered exten- 
sive subsequent deleptonization when passing through 
the neutrinosphere by electron captures on the shock- 
dissociated free protons. The low values of the lepton 
fraction (Y e u < 0.08) and the moderately high values of 
the entropy (s ~ 6) characteristic of this shocked ma- 
terial is the reason that ( -§^y- ) > 0, in accordance 

with Figure © and the discussion at the end of Section 
13.11 Finally, the material in this region tends to have a 
positive entropy gradient, a relic of the shock ramping 
up through this region. 

A positive gradient in s, positive cross response func- 
tions, and a positive derivative ( d\nY e ) can ' uno ^ er 

circumstances to be described, subject the region to an 
instability that we will dub "lepto-entropy fingers" (L- 
E F). This instability is thermal and lepton transport 
driven, as in the case of neutron fingers, but is very dif- 
ferent from the latter. Recall that the neutron finger 
mechanism, to operate in the presence of a positive (sta- 
bilizing) entropy gradient and a destabilizing lepton frac- 
tion gradient, requires that thermal transport be much 
more rapid than lepton transport. The argument given 
for the presence of the neutron finger instability in proto- 
super novae by the Liver more group, in fact, is based on 
the premise that thermal transport is rapid and lepton 
transport is negligibly slow. However, our results, which 
will be presented in detail in Sectional and which span an 
extensive grid of thermodynamic states and fluid element 
radii show that typically |TyJ ~ 10 x |E S |. From the def- 
initions of E s and Yy £ , given by equations (0), this in- 
dicates that lepton transport is considera bly more rapid 
than ther mal transport, a fact no ted by iBruenn et alJ 
(PHI) andlBruenn fc Dineval (119961). The re ason for this 
was pointed out in lBruenn fc Dineval ([1996ft and is that 
the energy, e v , transported by a neutrino is limited by the 
rapid rise in neutrino opacity with e v , whereas a v e or a 
9 e transports the same lepton number (viz., yi u = 1 and 
yi u = —1 for a v e and a v e respectively) independently 
of its energy. Furthermore, the v e and v e flows between a 
perturbed fluid element and the background can be op- 
positely directed, that is, additive in lepton number and 
subtractive in energy, depending on the temperature and 
the v e and v e chemical potential of the fluid element rel- 
ative to the background. 

Of critical importance in understanding the instabil- 
ity exhibited on the right-hand side of Figure [KJ the 
lepto-entropy finger instability, is the fact that both 

( din y £ ) > and the cross response function T s are 

large and positive in magnitude, typically as large or 
larger than |TyJ, and that the following inequalities are 
satisfied 

|Ti| » |Efc| (59) 

where z, j, and k can each denote either s or Yg. Thus, to 
operate in the presence of a positive (stabilizing) entropy 
gradient, lepto-entropy fingers requires lepton transport 
to be much more rapid than thermal transport, the oppo- 
site of the neutron finger requirement, and additionally 



requires that ^ d\nY g ) ^ ^ anc ^ ^ na ^ ^- s ^ e ^ ar S e and 
positive. 

The nature of the lepto-entropy finger instability can 
be illustrated by using inequalities to simplify equa- 
tions © - GU to give 



• ds 

dz 

= T S S + Ty,0- 



Y e VY t 



(60) 
(61) 



and 



(62) 



P\ ds Jp,Yi P\ dY tJ Pj s 

We have dropped the two small transport terms in equa- 
tion to get equation ([50)1 . and we have dropped 6y £ 
and the last term on the right-hand side of equation ([HI)) 
in comparison with the large and oppositely signed trans- 
port terms to get equation ([5T)l. Solving equation ([5T)l 
for 0y g gives 



M*) = 



(63) 



where we have written equation ([53)l explicitly for the 
case in which T s > and Ty £ < 0. Equation ([63)1 to- 
gether with equation ([55)) says essentially that the change 
of S given by equation ([55)) is quasistatic compared with 
the rate at which the large transport terms T S S and 
^Y^Yi can adjust #y?, and that 0y £ will therefore adjust 
to the instantaneous value of S as dictated by equation 
([55)) . 

Equation ([53)) is critical to understanding the lepto- 
entropy finger instability, and the lepto-entropy semicon- 
vective instability described in the next Section. While 
S arises in response to a displacement of the fluid ele- 
ment through a gradient in the background entropy, as 
dictated by equation ([55)) and by the identical equation 
([T5)) for the familiar Ledoux case, 0y £ arises in response 
to S , as dictated by equation ([53)) . rather than by a dis- 
placement of the fluid element through a gradient in the 
background lepton fraction, as dictated by equation {2J) 
for the Ledoux case. 

Plugging equation ([53)) into equation (|62|) gives 



<91n p 

ds 



<91n p 



■Yt 



(64) 



The derivative ( ) in equation ([5^)) is always < 0, 

but if the derivative (j^f^J > 0, which occurs in the 

outer region of the core where Yg is low as discussed above 
(e.g., Figure 01 and Figures E3 - EJ> , and if 



din p 
~8s~ 



P,Y £ 



<91n p 



> (65) 



then the fluid is unstable. For suppose, as an example, 
that ds/dz > 0. In the absence of transport this is a sta- 
bilizing gradient, for an outward displacement of a fluid 
element at constant entropy will result in 6 S < which 
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will make it more dense than the background and drive 
it back. Now, however, equation with equation (|65|) 
satisfied asserts that S < will lead to an outwards ac- 
celeration of the fluid element, i.e., the fluid is unstable. 
A more mathematical statement of the above can be ob- 
tained by using equation (|60|) for 6 S in the derivative of 
equation to get 



<91n p 

ds 



d\n p 



ds 



> 



dz 

(66) 

and if equation 



dlnp \ 

J pXe 



which has growing solutions if ^| 
is satisfied. 

How does this curious instability come about? A dis- 
placement of the fluid element through an entropy gra- 
dient results in a difference, # s , between the entropy of 
the fluid element and the background. The large value of 
T s means that the appearance of the entropy difference 
S will quickly induce via equation a lepton frac- 
tion difference, 0y i1 between the lepton fraction of the 
fluid element and the background. The 0y i so induced 
is dependent only on # s , has the same sign as # s , and 

is independent of the Yi gradient. Because 

and ^7^77^ are of opposite sign, if the perturbation S 

results in a restorative force, the induced perturbation 
6y £ will result in an anti-restorative force. If the induced 
perturbation 0y i is large enough relative to S1 the anti- 
restorative force wins out, and the perturbation grows, 
i.e., a "lepton finger" sustained by the entropy difference 
between the fluid element and the background will pene- 
trate into the background. Thus, a Ledoux stable region 
can be destabilized by this diffusive lepto-entropy finger 
instability. Crucial to the emergence of the lepto-entropy 
finger instability is the large positive value of T s , and we 
shall illustrate by an example in Section 0] how this large 
positive value arises. 

Despite the rather extreme assumptions that went into 
the derivation of inequality (|65|) (viz., |T S | >> uy n 
\T Yi \ » u Yn |E S | << u s , |£yj << u) s ), inequality 
(|65|) proves to be remarkably robust for predicting the 
lepto-entropy finger instability in proto-supernovae. In 
Section |S| we test this criterion for lepto-entropy fingers 
for one of the models we analyze by reducing T s to vio- 
late the inequality in equation by 5 percent for every 
core radius and fluid element size for which it was satis- 
fied, and find that the region previously found unstable 
to lepto-entropy fingers almost completely disappears. 

The above analysis of the lepto-entropy finger insta- 
bility breaks down when \ds/dz\ « \dYi/dz\. An out- 
ward displacement will then result in a very small mag- 
nitude of S and the last^ term in equation ([TU|) involv- 
ing the relatively large dYi/dz cannot now be neglected. 
Thus, there is a small sector about the vertical axis where 
\ds/dz\ « \dYi/dz\ which bounds the lepto-entropy fin- 
ger region in Figure ^| where the fluid is either stable or 
convective. 

3.4.2. Lepto-Entropy Semiconvection (L-E-Semi-C) 

Consider now the case for which ^| < (Ledoux desta- 
bilizing) but for which the conditions leading to equation 



(H2J are again satisfied (viz., |T S | >> uy n |TyJ >> uy n 
|E S | << u Sl |EyJ << u s ). Then equation (|66|) as it 
stands gives stable oscillatory solutions if equation JHSJ) is 
satisfied. However, the approximations leading to equa- 
tion (jnnj) are now inadequate, as the addition of even a 
small term in v, due for example to the samll but nonzero 
values of E s and Yly £ in equation ®, can, depending on 
its sign, turn the oscillatory solution of equation JSHI) 
into a damped oscillatory solution (stable) or a growing 
oscillatory solution (semiconvective) . Equation (JHEJ) was 
obtained by approximating equations (0 - (|TT|) by equa- 
tions (inOJ) - (|^ . To obtain a better approximation, we 
approximate equations (0 - {11} instead by 



Os — £,s#s 



= T,0, 



(67) 
(68) 



(69) 



*=-(¥) e -- g -(w) 

P \ ds J P ,Yi P 

Here we have assumed, as in our analysis of lepto-entopy 
fingers above, that equation ifTTijl is dominated by the 
transport terms and have dropped both the 9y e and the 
^j^v terms to get equation l|68jl . Unlike our analysis 
of lepto-entopy fingers, however, we have retained the 
transport terms in equation Q to get equation l|67|> . 
Solving equation lfo%|l gives as before 



1±. 



(70) 



where we have written the last expression in equation 
(|7n|) to explicitly account for the fact that T s > and 
Ty £ < 0. Here as before the large values of the response 
functions T s and Ty e in comparison with the terms we 
have dropped guarantees that the value of 0y £ adjusts 
immediately to the instantaneous value of S . Substitut- 
ing equation (|TU|) into equation (|69|) gives equation (|64|) 
again, i,e., 



<9ln p 

ds 



d\n p 



0s- (71) 



Taking the time derivative of equation {7TJ and using 
equation (|67|) for t9 s , we obtain in place of equation (|66|) 



din p 

ds 



P,Ye 



d\n p 



X [Tj s s + T l y £ 0y e 



ds 
dz 



d\np 
~ds~ 



P,Yt 



din p 



-Yt 



1 Y £ 

= (^s + Ey £ 

dlnp\ ( d\n p 



ds 
dz 



■Yt 



ds 



P,Y £ 



dYo 



■Y t 



ds^ 
dz 



v(72) 



16 



where we have used the equation (|T0|) in obtaining the 
second expression of equation (f72j) . and equation (|TT|) 
in obtaining the last. Equation (|72|l is a more accurate 
version of equation (|66j) and differs from the latter in that 
some of the terms multiplying v have been included. It 
will be used in place of equation (|66j) when the coefficient 
of v are needed to decide between growing or decaying 
oscillatory solutions. 

To continue our consideration of the case in which ^ < 
0, rather than examining equation (|72)l it is conceptually 
easier to examine its first integral 



Ey £ 



■Yt 



din p 

ds 



<91n p 



■§)<*■ 



We note first 
efficient of x 



from equation 
— xq is greater 



xo). 



(73) 



requires that ^ 



/ dlnp \ 



if the co- 
zero (i.e., 

> 0, which, of course, 



that 
than 



(dlnp\ 



< 0), then 



>0 since 

the solutions are oscillatory. The scenario here is analo- 
gous to that giving rise to lepto-entropy fingers discussed 
above. In fact, the same mechanism that destabilizes 
a fluid element in the presence of a stabilizing entropy 
gradient in the case above of lepto-entopy fingers, here 
tends to stabilize a fluid element in the presence of a 
destabilizing entropy gradient. Thus, an outward dis- 
placement of a fluid element through a negative entropy 
gradient results in the entropy of the fluid element be- 
coming greater than that of the background, that is, a 
S > develops. This by itself would lead to a posi- 
tive buoyancy force that would tend to force the fluid 
element farther outward. However, a S > induces a 
0y e > 0, in accordance with equation (|70j) . If 0y £ is large 

V ^ J P ,s 

the negative buoyancy arising from 0y g > together with 

a positive {^^f^j will overcome the positive buoyancy 

arising from S > 0, and the fluid element will be forced 
back. The fluid element will therefore oscillate. 

Whether the oscillatory motion of the fluid element is 
growing or decaying depends on the sign of the coefficient 
of v in equation (|73j) . If this coefficient is less than zero, 



enough, i.e., if 



(dlnp\ 

V ds J 



P,Yi 



> 0, then 



i.e., if E s 



< 0, the solutions of equation 



{73J) are decaying oscillatory and the fluid is stable. If 
> 0, 



on the other hand, the solution 

of equation (|T5|) is growing oscillatory and the fluid is 
unstable to semiconvection. 
To examine the sign of ^E s + Ey £ ) further, we 

note that the response function E s is always negative. 
If the cross response function T i y i were zero, the coeffi- 
cient of v in equation (|73|) would be negative, the solu- 
tions would be decaying oscillatory and the fluid would 
be stable. This is analogous to the Schwarzschild sta- 



bility described in Section with E s < 0. As pointed 
out there, an entropy perturbation, S , and the buoy- 
ancy force it induces will be diminished by the tendency 
of neutrino transport (modeled by a negative E s ) to ther- 
mally equilibrate the fluid element with the background. 
Here the situation is similar except that the buoyancy 
force arising from an entropy perturbation, # s , comes 
about through the 0y e that it induces. In either case the 
effect of a E s < is to reduce the buoyancy force with a 
phase such that a net work is extracted from the motion 
of the fluid element during the course of each oscillation. 

However, if Ey^ > and additionally satisfies the in- 
equality 



■Y t 



(74) 



then the coefficient of v in equation J73J) is > 0, the 
solutions are growing oscillatory, and fluid is semicon- 
vectively unstable. In this case, where both E s and Ey^ 
are nonzero, there are two competing effects governing 
the growth or decay of the amplitde of oscillation of the 
fluid element. These are embodied by the two terms in 
the coefficient of v in equation (|73|). The E s term, as de- 
scribed above, models the effect of an entropy difference 
S on the thermal transport between the fluid element 
and the background viz., for E s < the induced thermal 
transport drives a S tending to reduce the magnitude 
of 9 S , and with it the buoyancy force. The effect is that 
a net work is extracted from an oscillating fluid element 
during the course of each oscillation. The term involving 
Ey £ affects the motion more indirectly. It begins with the 
fact that an entropy difference, 9 S , induces a lepton dif- 
ference, 0Y e , in accordance with equation (|70j) . This lep- 
ton difference, being, in turn, another driver of thermal 
transport between the fluid element and the background, 
as modeled by equation ((67)). provides a contribution to 
6 S . This contribution has the same sign as S if Ey £ > 0, 
which is usually the case. If inequality (|74|) is satisfied, 
the contribution to 6 S arising from the term involving Ey^ 
is larger in magnitude than the contribution to S arising 
from the term E s , and the net effect is a 6 S having the 
same sign as S . With this, the oscillation amplitude of 
the fluid element grows. We thus have a case in which 
a lepton difference, 0y g1 induced via equation (JHH) by an 
entropy difference, 9 Sl does two things. It causes a grow- 
ing mode to become oscillatory (the coefficient of x — xq 
in equation {ZHJ)), and it causes this oscillatory solution 
to grow (the coefficient of v in equation (|73)0 . We refer 
to this instability as "lepto-entropy semiconvection," and 
denote it by "L-E-Semi-C" in the figures. 

The condition for lepto-entropy semiconvection ex- 
pressed by inequality ((71)) can be tested by reducing T s 
so that inequality (|73| is no longer satisfied. Doing so 
for the response functions shown in Figure El required 
a substantial reduction in T s before the lepto-entropy 
semiconvective instability disappeared. However, the 
magnitude of response functions in Figure El are not 
large enough for the velocity term in equation (P"U|) to 
be neglected in comparison with the diffusion terms. In- 
creasing all of the response functions by a factor of three 
(equivalent to reducing the size of the fluid element) gives 
the results shown in Figures El an d El Comparison of 
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these two figures shows that lepto-entropy semiconvec- 
tion is present for the case shown in Figure IT6l in which 
inequality (|74|) is satisfied, but not in Figure El m which 
the inequality is not satisfied. 
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Fig. 16. — Fluid stability as a fun ction of the Ins and lnYg 
gradients for the case shown in Figure 1151 but with each of the 
response functions increased by a factor of three. 
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Fig. 17. — Fluid stability as a function of the Ins and lnYg 
gradients for the values of the gravi tati onal acceleration and the 
response functi ons shown in Figure ITTTI but with T s reduced so 
that inequality 1741 is no longer satisfied. 



3.4.3. Lepto-Entropy Fingers and Lepto-Entropy 
Semiconvection in regimes where ^ §^py £ ) < 

Another important case for proto-supernova is the case 
in which the conditions leading to equation ((SI)) are again 
satisfied (viz., |T S | >> u)y t , |TyJ >> ujy n |E S | << uo s , 

l^yj << ^s)t but now (^^j < 0, which occurs in 
the cold inner core where Yg is high. Figure IT51 shows the 
stability /instability regions in the - d plane for 
a thermodynamic state representative of the colder and 
less highly deleptonized inner core where the above con- 
ditions are satisfied. The instabilities encountered here 



are similar to those exhibited in Figure and the dis- 
cussion here will therefore be applicable to the conditions 
shown in Figure ITU 

To examine this case, we note that the conditions 
stated above imply that the approximations leading to 
equation ((73} are satisfied here, and we rewrite this equa- 
tion as 











v = ^E s 








<91n 
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dY e 
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X ). 



(75) 



where we have explicitly exhibited the fact that now both 

(^r) Y < and (fff ) < °- E q uation G3 gives 
a stability criterion similar to the Schwarzschild stability 
criterion of equation (|55J). viz., stable oscillatory solu- 
tions if ^| > and unstable solutions if ^ < 0. In 
the latter case, we have convection if we are in a region 
which is Ledoux unstable (most of the left-hand side of 
Figure [jSt , an d lepto-entropy fingers in the region which 
is Ledoux stable (upper wedge on the left-hand side of 
Figure IT%|) . The explanation for the presence of lepto- 
entropy fingers here is similar in origin to the case dis- 
cussed in Section 13.4.11 Here a negative (destabilizing) 
gradient in s is stabilized, in the absence of transport, 
by a sufficiently large positive gradient in Yg depending 

on the relative magnitude of ( ) and ( ) 

With transport, however, the lepton fraction difference, 
0y l1 between the fluid element and the background no 
longer arises from the displacement of the fluid element 
in a background gradient of Y^, but in response to the 
entropy difference S in accordance with equation (|70[) . 
as in the case discussed in Section l5.4.1l Rather than sta- 
bilizing the fluid element, the induced 6y £ with the same 

( <91n 

sign as 6 S , together with the negative value of ( -^j 

further destabilizes the fluid element, and we have a case 
which we will again refer to as lepto-entropy fingers. 

The case in which ^| > is encountered frequently 
throughout much of the extent of the cold inner core 
of a proto-supernova. This positive entropy gradient 
arises during infall because the electron captures that 
take place in matter with increasing radius are more out 
of equilibrium due to a faster infall velocity of this mat- 
ter. (Recall that the inner core collapses approximately 
homologously with infall velocity proportional to radius.) 
This causes an increased entropy production with radius. 

A positive entropy gradient in equation {75} gives os- 
cillatory solutions. In this case the entropy gradient is 
stabilizing, i.e., an outward displacement of the fluid ele- 
ment gives rise to a negative entropy difference, # s , which 



/ p,s 



gives rise to a neg- 



together with a negative ( 

ative buoyancy forcing the fluid element inward. A neg- 
ative 0Y e arises in response to the negative 6 S , in ac- 
cordance with equation ((70)). and this together with the 

negative ( ) adds to the restoring force. 
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Fig. 18. — Fluid stability as a function of the Ins and lnYg 
gradients for the specified thermodynami c sta te for a case in which 
the response functions satisfy inequality 1741 . 
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Fig. 19. — Fluid stability as a function of the Ins and lnYg 
gra dients for the same thermodynamic state as shown in Figure 
1181 but for a c ase in which the response functions do not satisfy 
inequality 1741 . 



These oscillatory solutions are stable or semiconvec- 
tive depending on whether the sign of the coefficient of 
v in equation ((75)) is negative or positive, respectively. 
If the sign is positive, then the discussion at the end of 
Section 13.4.21 applies. In particular, an entropy differ- 
ence S drives a thermal diffusion which tends to reduce 
the magnitude of S , damping the oscillation of the fluid 
element. But the 6y e that arises in response to # s , in 
accordance with equation ((70)). feeds back by driving a 
thermal diffusion tending to increase the magnitude of 
S . If this latter thermal diffusion is large enough (i.e., if 
inequality ((71)) is satisfied) it dominates and causes the 
oscillation amplitude of the fluid element to grow. We 
then get a semiconvection driven by a dy t induced by a 
6 S , and we will refer to this case again as lepto-entropy 
semiconvection. 

Referring back to Figure ^1 the values of the response 
functions in this figure were chosen to satisfy inequality 
{73J) but to otherwise be representative. Note that the 



region on the right (where || > 0) is unstable to lepto- 
entropy semiconvection, in accordance with our discus- 
sion discussed above. As a test of the condition for lepto- 
entropy semiconvection as prescribed by inequality ((71)) . 

FigureHHIshows the - ^jp- plane for the same ther- 
modynamic state and response function values as Figure 
El except that the value of the response function T s has 
been decreased from 7 x 10 4 to 5 x 10 4 so that inequal- 
ity O is now violated. The region in the - ^fp- 
plane unstable to lepto-entropy semiconvection in Fig- 
ure El is now stable, validating the prescription given by 
inequality ((71)) for this case. 

4. COMPUTATION OF THE RESPONSE FUNCTIONS 

In the preceding Section we discussed the effect of dif- 
ferent sets of values of the response functions, S s ,Xly € , 
T s , and Ty £ , on the locations of the stability and in- 
stability regions in the ^j^- - dl ^ £ plane, and on the 
modes of instability. In this Section we discuss how we 
obtain these response functions for a given thermody- 
namic state and fluid element size. Because the response 
functions depend on details of neutrino transport, and 
may be sensitive to these details, our approach is to com- 
pute the response functions by sophisiticated radiation- 
hydrodynamical equilibration simulations. The thermo- 
dynamic conditions (i.e., the values of the density p, tem- 
perature T, and lepton fraction Yg) of spherical fluid el- 
ement are specified, and the size (radius) of the fluid 
element is chosen. The fluid element is then represented 
numerically by 20 spherical mass shells of equal width, 
and the background by 60 additional mass shells of the 
same thickness surrounding the fluid element. Initially 
the fluid element and the background are at the same 
thermodynamic conditions. The fluid element is then 
subjected to a perturbation in its thermodynamic condi- 
tions relative to the background, and the neutrino medi- 
ated re-equilibration of the fluid element with the back- 
ground is then numerically computed. The radiation- 
hydrodynamical code used to perform these equilibra- 
tion simulations is the same as that used to model core 
collapse and will be described in the next Section. 

Two procedures, "isothermal" and "adiabatic," for 
perturbing a fluid element relative to the background 
were tried. An isothermal perturbation proceeds as fol- 
lows. Prior to performing the perturbation, the thermo- 
dynamic conditions of the fluid are frozen, neutrino in- 
teractions with the matter, but not transport, are turned 
on, and neutrinos are allowed to equilibrate locally with 
the fluid element and the background. The result is a 
uniform neutrino sea of all flavors in equilibrium with 
both the fluid element and the background. To com- 
plete the isothermal perturbation, the neutrino distribu- 
tions are kept constant and the temperature or electron 
fraction or both of the fluid element are perturbed rel- 
ative to the background, with the density of the fluid 
element adjusted so that pressure equilibrium between 
the fluid element and the background is maintained. In 
summary, an isothermal perturbation is a perturbation 
of the thermodynamic state of the fluid element relative 
to both the background and a uniform neutrino "sea." 
To complete the simulation, neutrino transport is turned 
on. As the neutrinos re-equilibrate the fluid element with 
the background the radius, and therefore density, of the 
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fluid is continually adjusted so that pressure equilibrium 
is maintained at all times. The results of this simulation 
are the time evolutions of the quantities (s) — s = S and 
(Yi) —Yi = 0y 1 - ) where (s) and (Yi) are the entropy and 
lepton fraction mass averaged over the fluid element, s 
and Yi are the entropy and lepton fraction of the back- 
ground, and S and e Ye are defined at the beginning of 
Section □ 

An adiabatic perturbation begins with a perturbation 
of the temperature or electron fraction or both of the 
fluid element relative to the background. To complete 
the adiabatic perturbation, neutrino interactions with 
the matter, but not transport, are turned on, and neu- 
trinos are allowed to equilibrate locally with the fluid 
element and the background. The result is that both the 
fluid element and the neutrinos within it are perturbed 
with respect to the background and the background neu- 
trinos. Thus, an adiabatic perturbation is a perturbation 
of both the matter and neutrinos such that local thermo- 
dynamic equilibrium between the matter and neutrinos is 
maintained. Once the neutrinos have achieved local equi- 
libration with the matter, neutrino transport is turned 
on and the re-equilibration of the fluid element with the 
background is followed. Pressure equilibrium between 
the fluid element and the background is maintained con- 
tinuously by adjusting the radius (and therefore density) 
of the fluid element. As before, the results of this simula- 
tion are thejtime evolutions of the quantities (s) — s = S 
and (Y £ )-Y £ = 0Y t . 

From a given simulation, time derivatives of the quan- 
tities S and 0y i are computed by taking finite differences, 



i.e., 



0v £ = 



8 (A 8 t) - 8 (O) 
A s t 

e Yi (A Yi t) - o Yi (o) 



(76) 



where we have defined A s t and A Yi t to be the times 
required for S and 6 Ye to be reduced to 1/e of their orig- 
inal perturbed values, respectively. Experiments showed 
that the choice of 1/e as a measure of the approach to 
equilibration in the determination of S and Yi was not 
critical, although the use of shorter time intervals gave 
slightly higher values for S and 6 Ye , as would be ex- 
pected. 

To determine values of the response functions for a 
given thermodynamic state and fluid element radius, we 
performed two equilibration simulations beginning with 
different linearly independent sets of values of S and Ye 
for the perturbation of the fluid element with respect 
to the background. Call these two sets of perturbations 

Oyj and oi 2 \Oy^- We then took a linear combina- 
tion of the two perturbations, choosing the coefficients a 
and b such that 



0s = ae a 1 ^+b0 a 2 ^O,0 Yi 



0, (77) 



and another linear combination, choosing c and d such 
that 



sflW+dfl?) =0, 



c0 



(i) 

Y e 



d9^ + 0. (78) 



From these linear combinations we compute the response 
functions. For example, 
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CU Y e 



(80) 



(Naturally it would have been possible to perturb 6 S 
at e Yl = 0, and vice versa, and this would have made it 
unnecessary to chose linear combinations of the results. 
However, the procedure we adopted proved to be simpler 
to implement.) 
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Fig. 20. — Temperature as a function of radius at selected times 
for an equilibration simulation with improved neutrino rates. The 
radius of the fluid element is 1 km, and the thermodynamic state 
of the background and the selected times are indicated. 

As we discussed above in Section a large positive 
value of T s is critical to the emergence of the lepto- 
entropy finger instability. To illustrate how this comes 
about with a specific example, we show in Figures I20l-l25l 
some of the details of an equilibration simulation for an 
isothermal perturbation. The simulation was performed 
with improved neutrino rates. For ease of illustration, 
the state of the background that we chose, p = 3xl0 12 g 
cm -3 , s = 7, Yi, is a thermodynamic state for which 
the v e (and therefore v e ) chemical potential vanishes. 
The equilibrium v e and v e number densities in the back- 
ground are therefore equal to each other. Because of 
this, it was possible to perturb the temperature of the 
fluid element and adjust its density to keep the pressure 
the same as the background while maintaining Yg at its 
unperturbed value. (Different values of v e and v e num- 
ber densities would have resulted in a perturbed value of 
Yi in the fluid element after its density was adjusted to 
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Fig. 21. — Same as Figure 1201 but showing the electron fraction 
as a function of radius at selected times. 



Fig. 23. — Same as Figure 1201 but showing the entropy as a 
function of radius at selected times. 



maintain pressure equilibrium, even if its value of Y e was 
kept equal to that of the background.) While the choice 
of any other thermodynamic state would have illustrated 
our points, a choice which makes it possible perturb the 
temperature (and entropy) of the fluid element, main- 
tain pressure equilibrium, and keep the value of Yg at 
its unperturbed value provides the simplest and cleanest 
illustration of how a lepton fraction perturbation arises 
in response to an entropy perturbation, which is the ba- 
sis of both the lepto-entropy fingers and lepto-entropy 
semiconvection instabilities. 
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Fig. 22. — Same as Figure l20l but showing the entropy as a 
function of radius at selected times. 



Figures EH and El show the matter temperature and 
electron fraction profiles at selected times. The fluid 
element was perturbed relative to the background by 
decreasing its temperature, leaving the electron frac- 
tion unchanged, and adjusting the density to restore 
pressure equilibrium with the background. This initial 
perturbation is shown in the figures by the tempera- 
ture and electron fraction profiles at the smallest time, 
t = 1.44 x 10 -10 . The above perturbation in T results 



-10 



in the entropy perturbation shown in Figure [221 by the 
profile at t = 1.44 x 10 -10 . Because of our choice of 
a thermodynamic state with zero v e chemical potential, 
there is no net perturbation in the lepton fraction, as 
shown in Figure 1251 bv the profile at t = 1.44 x 10 -10 . 
As explained earlier in this Section, an isothermal per- 
turbation is effected by first allowing the neutrinos to 
equilibrate locally with the fluid element and the back- 
ground, both at the same thermodynamic state. The 
matter, but not the neutrinos, of the fluid is then per- 
turbed. Thus, in particular, the electron neutrino and 
antineutrino number densities, n Ve and n^ e , are initially 
uniform when the matter is perturbed, as shown by the 
n Ve and n Pe profiles in Figures El and [23 respectively, at 
t = 1.44 x 10" " 
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Fig. 24. — Same as Figure l20l but showing the v e number density 
as a function of radius at selected times. 



What happens almost immediately after the perturba- 
tion is that the ^ e 's and zVs in the perturbed fluid ele- 
ment equilibrate locally with the matter, essentially es- 
tablishing, thereby, conditions equivalent to an adiabatic 
perturbation. This is evident by comparing the v e and 
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Fig. 25. — Same as Figure l^Ol but showing the u e number density 
as a function of radius at selected times. 



v e profiles at times 1.7 x 10 _T s and 1.06 x 10 -6 s in Fig- 
ures and with the corresponding profiles in Figures 
EE! and El . These latter figures show the time evolution 
of local the neutrino-matter equilibration with transport 
turned off so that only local equilibration can occur. It is 
seen that the profiles with transport turned on (Figures 
El and ESI) and the corresponding profiles with transport 
turned off (Figures El and 127)1 follow each other closely, 
and that local equilibration is achieved in about 1 /i s. 
After this time the effects of transport begin to be seen, 
and the subsequent evolution is the same whether the 
initial perturbation is isothermal or adiabatic. This is 
why there is very little difference in the values of the re- 
sponse functions whether the perturbation is isothermal 
or adiabatic. 

The local equilibration is very asymmetric, with n Pe 
decreasing by ~ 12.5 % versus a ~ 0.36 % decrease in 
n Ue , and this asymmetry is responsible for much of what 
follows. The asymmetry is due to the fact that the equi- 
librated values of n Ve and np e depend on electron num- 
ber density, n~, and positron number density, n e +, re- 
spectively, through the reactions e~ + p ^ n + v e and 
e + + n ^ p + V e . The positrons are thermally produced 
as electron-positron pairs, so the value of n e + is very 
sensitive to temperature and drops significantly with the 
temperature perturbation, whereas n~ consists largely of 
fixed electrons (for the entropies of interest here), and is 
much less sensitive to temperature. Thus, while the per- 
turbation of the fluid element reduces both n~ and n e + 
equally, the relative reduction in n e + is much greater. 
The result is that the local equilibration of the z/ e 's and 
zVs with the matter entails more v e absorption on pro- 
tons (Figure El profiles at 1.7 x 1(T 7 s and 1.06 x 1(T 6 
s) than v e absorption on neutrons (Figure EJ profiles at 
1.7 x 10 -7 s and 1.06 x 10 -6 s), and this depresses Y e 
(Figure El profiles at 1.7 x 1(T 7 s and 1.06 x 1(T 6 s). 

After local equilibration, rapid lepton transport be- 
tween the fluid element and the background adjusts the 
n Ue and n Pe distributions toward values for which there 
is a quasi net lepton flow equilibrium (i.e., a small ratio 
of net lepton flow to the flow in each direction) of z/ e 's 
and VqS. Since the local matter- neutrino equilibration 



reduced the value of rip e much more than that of n Ve , the 
difference between these two quantities drives a net lep- 
ton flow out of the fluid element, raising the value of n Pe 
and reducing that of n Ve (FiguresElandEI respectively, 
profiles at 6.55 x 10 -6 s, 5.94 x 10 -5 s, and 1.16 x 10 -4 s) 
until the establishment of a quasi net lepton flow equilib- 
rium. The net flow of leptons out of the fluid element as 
a quasi net lepton flow is established reduces the value of 
Yg in the fluid element ( Figure l23l profiles at 6.55 x 10 -6 
s, 5.94 x 10" 5 s, and 1.16 x 10" 4 s). At this point, then, 
there is a rapid reduction in the value of 0y £ in response 
to the initial negative perturbation, # s , of the fluid ele- 
ment. This is the origin of the large positive value of T s , 
in accordance with its definition by equation El which is 
so crucial for the appearance of lepto-entropy fingers and 
lepto-entropy semiconvection. 

After the establishment of a quasi net lepton flow equi- 
librium, the differences, An Ue and Anp e , between the val- 
ues of n Ue and n Pe , respectively, between the fluid element 
and the background adjust as energy diffuses between 
the fluid element and the background, and slowly ap- 
proach zero as the temperature of the fluid element equi- 
librates with that of the background on a relatively long 
timescale (FigureEl profiles at 1.16 x 10 -4 s, 2.00 x 10 -4 
s, 5.58 x 10" 4 s, and 1.01 x 10" 3 s). 

In summary, an isothermal perturbation of the fluid el- 
ement is followed by a local equilibration of matter and 
neutrinos on the shortest timescale, producing conditions 
equivalent to that of an adiabatic perturbation. Rapid 
lepton transport then adjusts the values of n Ue and n Pe 
on the next shortest time scale, establishing a net lepton 
flow equilibrium and producing thereby a lepton fraction 
difference, 6y e , in response to the initial entropy differ- 
ence, 6 S . The slower thermal diffusion rate equilibrates 
the temperature of the fluid with the background on the 
longest time scale, while the lepton flow continually ad- 
justs the values of n Ue and n Pe to maintain a quasi net 
lepton flow equilibrium between the fluid element and 
the background. 
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Fig. 26. — Same as Figure EJbut with zone to zone transport 
turned off. 
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ideal gas (with excited states), and a nine-species nu- 
clear reaction network was used to follow the nuclear 
transmutations. We will refer to this equation of state 
as the "non-NSE" equation of state. A zone not in 
NSE were flashed to NSE if its temperature reac hed 0.44 
MeV ((Thielemann. Nomoto. fc Hashimoto1ll996ft . In or- 
der to compare our results more closely with those of 
the Livermore group, we also ran some of the simula- 
tions, de signated by " W," using the Livermore equation 
of state ( Wilson 2003). This equation of state was pro- 
vided in tabular form. This W equation of state was 
joined smoothly to the non-NSE equation of state, and 
the latter was used when nuclei were not in NSE. Tabled 
summarizes the core collapse simulations that were per- 
formed to produce the core configurations that will be 
analyzed in the next Section. 



Fig. 27. — Same as Figure 1^1 but with zone to zone transport 
turned off.. 



As representative of the cores of supernova progeni- 
tors, we have chosen to examine at 100 ms and 200 
ms after bounce the cores of a 15 M (S15s7b) and 
a 25 Mf?) (S25s7b) s tellar progenitor. These were pro- 
vided by IWooslevl (^995) and evolved from the main 
sequence to the on s et of core collapse as described in 
iWooslev fc Weaverl ((19951) . Precollapse model S15s7b 
has an iron core mass of 1.278 M , representative of 
a small iron core mass progenitor, while model precol- 
lapse model S25s7b has an iron core mass of 1.770 Mq, 
representative of a large iron core mass progenitor. Each 
progenitor was evolved through core collapse, bounce, 
and the shock rehea ting epoch by the code de scribed in 
iBruennl (|l985t ) and iBruenn fc HaxtQ"nl (Il99lf). and ex- 
tens ively revised a s described in iBruenn et alJ ((200 it ) 
and IBruennl ((20031) . The code is fully general relativis- 
tic. Neutrino transport is multigroup, three flavor and 
uses a flux-limited diffusion approximation. Twenty en- 
ergy zones were used in these simulations spanning in 
geometric progression the neutrino energy range from 
4 MeV to 400 MeV. T he n eutrino microphys i cs de- 
scribed in IBruennl ((19851) and IBruenn fc Haxtonl ((19911) 
was used for the simulations, as well as some re- 
vised and new neutrino microphysics described below. 
Also used were the ion screening corrections given by 
iHorowitzl (1 19961) and implement ed as described in detail 
in IBruenn MezzacapjSl l|1997D- 

In some of the simulations, designated by "LS," the 
Lattimer-Swesty equation of state (Lattimer & Swesty 
1991) was used when the following three conditions were 
satisfied locally: (1) n B > 10~ 8 fm~ 3 (p > 1.67 x 10 7 
g cm -3 ), where tib is the number density of nucleons 
(free and bound) per cubic Fermi, (2) T > 0.05 MeV, and 
(3) the matter was assumed to be in nuclear statistical 
equili brium (NSE). The Cooperstein-BCK equation of 
state ((Coopersteinlll985l:lBaron. Cooperstein. fc Kahanal 
119851) was used when the second and third condition 
was satisfied, but not the first. For nuclei not in NSE, 
which comprised the silicon layer, oxygen layer, and 
other exterior layers until they encountered the shock, 
seven alpha particle nuclei (plus neutrons and protons) 
were tracked explicitly. The nuclei were treated as an 



TABLE 2 
Core Collapse Simulations 



Simulation 


Progenitor 


Equation of State 


M15_LS 
M15_W 
M25_LS 
M25_W 


S15s7b 
S15s7b 
S25s7b 
S25s7b 


Lattimer-Swesty 
Livermore 

Lattimer-Swesty 
Livermore 



For each of the models included in Table [21 the the sta- 
bility of the fluid from the core center to the region just 
beneath the neutrinosphere was analyzed for two time 
slices, 100 ms and 200 ms after core bounce. The re- 
sponse functions, £ s ,£y £ , T s , and Ty £ , were obtained 
by performing equilibration simulations as described in 
Section El by the same code described above for the core 
collapse simulation. For most of the equilibration simula- 
tions, the "sta ndar d" neutrino microphysics d escribed in 
IBruennl (|l985l) and IBruenn fc Haxtonl (|l99ll ) was used, 
i.e., the same neutrino microphysics as used in the core 
collapse simulations. We observe, however, that neu- 
trino physics has advanced beyond this standard. Re- 
cently, analytic expressions have been derived for neu- 
trino absorption, emission, and elastic scattering on nu- 
cleons which take into account the nucle on recoil and 
block ing; factors for arbitrary degene racy (|Reddv et alJ 
119981: iBurrows fc Sawverlll9M 119991) . These new rates 
become smaller at high densities than the correspond- 
ing rates without nucleon blocking, and the recoil ef- 
fects call into question the assumption that the scat- 
tering rates are isoenergetic. The latter will predomi- 
nantly affect the thermalization of the z/ x 's (by V x 's" 
we refer to v^s and ^ r 's and their antiparticles), as 
these particles lack the energy exchange channels of 
the charged-current processes on baryons which dom- 
inate the thermalizatio n of the zVs and zVs . Ther- 
malization simulations ((Thompson et alJ 12000). for ex- 
ample, show that ^ x -nucleon scattering is an impor- 
tant thermalization process for z/ x 's with energies ex- 
ceeding 15 MeV. Other neutrino processes that may 
be important for equilibrating a fluid element with 
the background are neutrino-nucleon bremsstrahlung 
and the related neutrino-nucleon inelastic scattering 
(lHannestad fc RaffeltJIT993 iThomnson e t, al 1 12000ft . 

We have performed equilibration simula tions for sev- 
eral of the models with the improved (|Reddv et alJ 
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Il998|) absorption, emission, and elastic scattering on nu- 
cleons, and with the inclusion of bremsstrahlung and 
neutrino-nucleon ine l astic scattering as formulated by 
lHannestad fc Raifeltl (|1998T ). Models analyzed for sta- 
bility using this improved neutrino microphysics are de- 
noted by "Impr" , while those using the standard neutrino 
microphysics are denoted by "Stnd." Table |3| lists the 
models analyzed for stability. We did not deem it neces- 
sary to include models obtained from core collapse simu- 
lations employing the improved neutrino physics because 
they differ very little below the neutrinosphere from those 
computed with the standard neutrino physics up to 200 
ms after bounce. The improved neutrino physics mainly 
affects the luminosities and RMS energies of the z/ x 's near 
and above the neutrinosphere. Furthermore, we wanted 
to perform the core collapse simulations with the Liver- 
more equation of state and comparison models using the 
Lattimer-Swesty equation of state with neutrino micro- 
physics as similar to that used by the Livermore group 
as possible. 

6. THE STABILITY OF COLLAPSED CORES BELOW THE 
NEUTRINOSPHERE 

Before analyzing the models listed in TableOlfor stabil- 
ity using the method described in the preceding Sections, 
let us first compare the core profiles produced by simu- 
lations using the Livermore equation of state with those 
produced by the Lattimer-Swesty equation of state, and 
ascertain the Ledoux unstable regions and the regions 
unstable to neutr on fingers by the crite ria given by the 
Livermore group ((Wilson fc Mavlell988D . Recall that ac- 
cording to the Livermore group's criteria a fluid is neu- 
tron finger unstable if a fluid element displaced outwards 
at constant composition but in temperature and pres- 
sure equilibrium with the background ends up after the 
displacement with a density less than that of the back- 
ground. The profiles in entropy and lepton fraction of 
models M15_W and M15JLS are shown at 100 ms and 
200 ms after bounce in Figures EEI and OH respectively. 
Likewise, the profiles of models M25_W and M25JLS are 
shown at 100 ms and 200 ms after bounce in Figures EH 
andlSTl respectively. 

Figures 1251 - I3D show that the core collapse simulations 
using the Livermore equation of state produce profiles in 
entropy and lepton fraction at 100 and 200 ms similar 
in shape to the simulations using the Lattimer-Swesty 
equation of state. The main differences are the slightly 
higher overall values of the entropy and lepton fraction 
given by the Livermore equation of state. These differ- 
ences arise mainly during infall because of the smaller 
free proton fraction given by the Livermore equation of 
state in comparison with the Lattimer-Swesty equation 
of state for a given thermodynamic state. A smaller free 
proton fraction during infall leads to a reduced rate of 
electron capture, and this results in a higher lepton frac- 
tion at trapping. The electron captures that occur after 
trapping as the z/ e 's equilibrate with the matter produces 
a greater increase of entropy, as the smaller free proton 
fraction and reduced electron capture rate given by the 
Livermore equation of state results in a greater disequi- 
librium between matter and z/ e 's at trapping. 

All models exhibit a region of Ledoux instability, in- 
dicated by the checkerboard patterns in the horizontal 
stripes, in the vicinity of 20 km from the core center at 
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Fig. 28. — Entropy and lepton fraction profiles of models Ml 5_W 
and M15_LS (as described in Table 100 ms after bounce. Also 
shown are the Ledoux unstable regions (checkerboard hatching) 
and the neutron finger unstable regions (diagonal pattern) as given 
by the Livermore group's criteria. 
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Fig. 29. Same as Figure 1231 but for models M25_W and 
M25_LS. 



100 ms, and models M15_W and M25_W exhibit an addi- 
tional region of Ledoux instability in the vicinity of 8 km 
from the core center. At 200 ms the inward advection of 
matter has moved these Ledoux unstable regions slightly 
inward, and additional (very mild) Ledoux unstable re- 
gions have developed in models M15_LS and M25JLS at 
small radii. In most cases these Ledoux unstable regions 
are driven by negative entropy gradients. The location 
and even existence of these Ledoux unstable regions must 
not be taken very seriously, however, as convection, not 
included in these simulations, will rapidly smooth the 
entropy and lepton fraction gradients driving the con- 
vection. 

Also shown in these figures by the slanted line pattern 
in the horizontal stripes are the regions that are neu- 
tron finger unstable by the Livermore group's criteria. 
The neutron finger instability considerably extends the 
regions of fluid instability to either side of the Ledoux 
unstable regions. In particular, the region of fluid insta- 
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TABLE 3 
Models Analyzed for Stability 



Model 


Progenitor 


Time from Bounce 


EOS 


z/-Microphysics 


M15_LS_100_Stnd 


S15s7b 


100 


ms 


Lattimer-Swesty 


Standard 


M15_LS_100_Impr 


S15s7b 


100 


ms 


Lattimer-Swesty 


Improved 


M15_W_100_Stnd 


S15s7b 


100 


ms 


Livermore 


Standard 


M15_LS_200_Stnd 


S15s7b 


200 


ms 


Lattimer-Swesty 


Standard 


M15_W_200_Stnd 


S15s7b 


200 


ms 


Livermore 


Standard 


M25_LS_100_Stnd 


S25s7b 


100 


ms 


Lattimer-Swesty 


Standard 


M25_LS_100_Impr 


S25s7b 


100 


ms 


Lattimer-Swesty 


Improved 


M25_W_100_Stnd 


S25s7b 


100 


ms 


Livermore 


Standard 
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Fig. 30. — Same as Figure |28|but at 200 ms after bounce. 
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Fig. 31. — Same as Figure |2H but at 200 ms after bounce. 



bility extends nearly up to the neutrinospheres, which 
are located near the outer parts of the troughs in the 
lepton fraction profiles. This illustrates the claim by the 
Livermore group that fluid motions resulting from the 
neutron finger instability will advect lepton rich matter 



to the neutrinosphere from regions below, and thereby 
enhance the u e luminosity. 

Let us now consider the analysis of the models listed 
in Table 03 for stability using the methodology developed 
in the preceding Sections. To perform this analysis, each 
model was considered in turn and a double loop over the 
model was performed. The outer loop was over the spher- 
ical mass shells between the core center and the neu- 
trino sphere (which comprised between 55 and 70 mass 
shells, depending on the model). For each mass shell 
the thermodynamic state and the gravitational accelera- 
tion were noted, and an inner loop was performed over 
the fluid element size. For this inner loop the radius of 
the fluid element was set initially to the v e flux- averaged 
mean free path for the mass shell under consideration 
and then increased successively by a factor of a/10 as the 
loop was executed until the radius of the fluid element 
exceeded one-half the distance of the mass shell from the 
center of the core. When the latter occurred, the loop 
was completed by setting the fluid element equal to one- 
half the distance of the mass shell from the center of the 
core. Thus, fluid element radii were considered spanning 
the range (A^ e ) < i^biob < \R, where (A^ e ) is the flux- 
averaged v e mean free path, i^biob is the radius of the 
fluid element, and R is the distance from the core center 
to the mass shell in question. Executing this inner loop 
for each mass shell typically involved 6 to 12 fluid ele- 
ment radii. For each mass shell and fluid element radius 
two liinearly independent equilibration simulations were 
performed as described in Section^ The results of this 
pair of equilibration simulations were the four response 
functions, E s , Ey £ , T s , and Ty^. Having computed these 
response functions, the gradients of Ins and In at the 
mass shell in question were noted, and the roots of equa- 
tion (|53|) were obtained. 

From these roots the stability of the fluid element of the 
given radius and at the given mass shell was ascertained 
as follows. If all the roots were real and negative, or if 
a real root was negative and the real part of a complex 
conjugate pair of roots was also negative, then the fluid 
was deemed stable for the given fluid element radius. If 
one or more of the real roots was positive, or if a real 
root and/or the real part of a complex conjugate pair of 
roots was positive, then the fluid was deemed unstable 
with a growth rate given by the largest of the positive 
elements, i.e., the largest of the real roots or the largest 
of the real root and the real part of a complex conju- 
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gate pair, as the case may be. If the largest positive 
element was the real part of a complex conjugate pair, 
and lepton transport dominated thermal transport (i.e., 
|TyJ > |E S |. which was always the case), the fluid was 
deemed unstable to lepto-entropy semiconvection for the 
given fluid element radius. If the largest positive element 
was a real root, the fluid would otherwise be stable in the 
absence of thermal or lepton transport (i.e., Ledoux sta- 
ble), and lepton transport dominated thermal transport 
(again this was always the case), the fluid was deemed 
unstable to lepto-entropy fingers for the given fluid ele- 
ment radius. If the largest positive element was a real 
root, and the fluid would otherwise be unstable in the 
absence of thermal or lepton transport (i.e., Ledoux un- 
stable), the fluid was deemed unstable to convection for 
the given fluid element radius. 



M1 5_LS_1 00_Stnd Lepto-Entropy 

Semiconvection 





M15 LS 100 Stnd 



Convection 




Fig. 33. Model M15_LS_100_Stnd: Growth rate as a 
function of core radius and fluid element size for Convection. 

Figures 03 - E3 show the results of our analysis of the 
stability of model M15_LS_100_Stnd. Figure g3 gives 
the growth rate as a function of core radius and fluid el- 
ement size for Ledoux convection, i.e., the growth rates 



Fig. 35. Model M15_LS_100_Stnd: Growth rate as a 
function of core radius and fluid element size for lepto-entropy 
fingers. 



of the fluid instability that would aeise in the absence of 
thermal and lepton transport. Except for an extremely 
narrow region near the core center, only the region be- 
tween R ~ 17 km and R ~ 22 km is Ledoux unstable, 
where R is the core radius, as is also evident from Figure 
(EE}- 

Figure 0S1 gives the growth rates for the convectively 
unstable region, i.e., the region that is unstable in the 
presence of thermal and lepton transport that is also 
Ledoux unstable. The unstable convectively unstable re- 
gion is essentially the same as the Ledoux unstable re- 
gion, except that thermal and lepton transport slightly 
reduces the growth rates. This is because the convec- 
tively unstable region is buoyancy driven, and the ten- 
dency of thermal and lepton transport to equilibrate the 
fluid element with the background reduces the buoyancy 
that drives this convection. 

Most interesting are Figures and ISol which show re- 
gions that are unstable solely by virtue of thermal and 
lepton transport. Figure shows that the lepto-entropy 
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semiconvectively unstable regions are located toward the 
inner part of the core, i.e., between R ~ 8 km and R ~ 16 
km. Referring to the discussion in Section R .4 . HI and Fig- 
ure ^1 instability to semiconvection is anticipated in the 
cold and relatively lepton rich inner core when the en- 
tropy gradient is positive, as it is for 8 > R > 16 according 
to Figure E3 The lepto-entropy semiconvection will be 
dominated by small scales, as the growth rate for this 
instability increases with decreasing fluid element size, 
as is evident from Figure 021 Figure shows the re- 
gion unstable to lepto-entropy fingers. Significantly, it 
extends from the Ledoux unstable region to the vicin- 
ity of the neutrinosphere, and could be a very important 
source of lepton transport. The growth rates peak for 
fluid element radii ~ 1/20 x R near the neutrinosphere 
increasing in comparison with R to ~ 1/4 x R deeper in 
the core. 



M15_LS_100_Stnd 



Lepto-Entropy Fingers 
Y = < 0.95 Critical Value 




Fig. 36. Model M15_LS_100_Stnd: Growth rate as a function 
of core radius and fluid element size for lepto-entropy fingers when 
T s is limited to 95 pe rcen t of its critical value, as given by the 
inequality in equation 1651 . 



As a test of the applicability of inequality in pre- 
dicting the lepto-entropy finger instability in the col- 
lapsed cores of supernova progenitors, Figure OH shows 
that when the cross response function T s is restricted 
to be no larger than 95 percent of the critical value, 
as given by the inequality in equation ((65)). the lepto- 
entropy finger instability disappears except for an in- 
significant region near the core center. Figure OU pro- 
vides a global view of the significance of inequality 
for a particular core collapse model. Figures ?? and| 
provide an expanded view of the critical role of inequal- 
ity (|65[) at a given thermodynamic state. These latter 

two figures show the — plane for the ther- 

modynamic state (p = 3.4 x 10 12 g cm -3 , T = 10.3 
MeV, Y e = 0.1165), which pertains to one of the zones 
of Model M15_LS_100_Stnd. The ratio of fluid element 
size to distance to core center (i.e., Rbiob/Radius) was 
chosen to be 0.047, for which the lepto-entropy growth 
rate was maxmimum. Shown is the convective regions 
(red), lepto-entropy semiconvective region (green), lepto- 
entropy finger region (purple) and stable region (blue). 
Also shown are iso-growth rate contours. Figure 1551 
FigurelS^lshows the regions of convective instability for 
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Fig. 37. The 



d 1 " *i — d ^ s plane for the thermodynamic 
state = 3.4 x 10 12 g cm" 3 , T = 10.3 MeV, Y e = 0.1165), of one 
of the zones of Model M15_LS_100_Stnd, for which inequality J65j 
is satisfied, and Rbiob/Radius = 0.047. The color coding is red: 
convective, blue: stable, green: lepto-entropy semiconvective, and 
purple: lepto-entropy fingers. 
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Fig. 38. — The same as Figure ?? except that T s was r estri cted 
to be 95 percent of the value needed to satisfy inequality 1651 . 



model M15_W_100_Stnd, a model evolved with the Wil- 
son EOS. There is a convectively unstable region deep in 
the core (7 < R < 9 km), and another very narrow region 
near the core center. The convectively unstable region at 
7 < R < 9 km, not seen in model M15_LS_100_Stnd, corre- 
sponds to the negative gradients in both s and Yi, which 
can be seen in Figure (|28j). and which causes this re- 
gion to be Ledoux unstable. Figures H01 and HT1 show the 
regions of the lepto-entropy semiconvective and lepto- 
entropy finger instabilities, respectively. These regions 
are located in the core very similarly to the locations 
of the corresponding regions in model M15_LS_100_Stnd, 
except that the lepto-entropy growth rate is down by a 
factor of 2 or 3. 

Figures 02| and QHI show the regions of lepto-entropy 
semiconvective and lepto-entropy finger instabilities, re- 
spectively, for model M15_LS_100Jmprv, a model whose 
regions and stability and instability were analyzed on 
the basis of "improved" neutrino microphysics. Lepto- 
entropy semiconvective is dominated at small scales, as 
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Fig. 39. Model M15_W_100_Stnd: Growth rate as a function 
of core radius and fluid element size for Convection. 



Fig. 42. — Model M15_LS_100_Imprv: Growth rate as a function 
of core radius and fluid element size for semiconvection computed 
with "improved" neutrino rates. 
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Fig. 40. Model M15_W_100_Stnd: Growth rate as a function 
of core radius and fluid element size for semiconvection. 



Fig. 43. — Model M15_LS_100_Imprv: Growth rate as a func- 
tion of core radius and fluid element size for lepto-entropy fingers 
computed with "improved" neutrino rates.. 
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Fig. 41. — Model M15_W_100_Stnd: Growth rate as a function 
of core radius and fluid element size for lepto-entropy fingers. 



in the case of the standard neutrino microphysics (Fig- 
ure but extends over a smaller region of R, viz., 
S^R^12 km. Lepto-entropy fingers extends over a 
larger range of R, from R = 15 km up to the vicinity of 
the neutrinosphere, instead of from 23 km in the case of 
the standard neutrino microphysics. The lepto-entropy 
finger growth rates are about the same in both analyses, 
being about several hundred s -1 throughout most of the 
unstable region. 

Figures ^] - 03 show the various instability re- 
gions for model M25_LS_100_Stnd and Figures SHI - 
E3 do the same for M25_W_100_Stnd (we omit the 
"Ledoux Convection" graph here because it is very sim- 
ilar to the "Convection" graph). Figures and 
show the regions of lepto-entropy semiconvective and 
lepto-entropy finger instabilities, respectively, for model 
M25_LS_100 Jmprv (we omit the graphs of Ledoux Con- 
vection and Convection because they are very simi- 
lar to the results for model M25_LS_100_Stnd. The 
graphs for models M25_LS_100_Stnd, M25_W_100_Stnd, 



28 



M25 LS 100 Stnd 



Ledoux Convection 








M25_LS_100_Stnd 



Lepto-Entropy Fingers 




rr 



Fig. 44. — Model M25_LS_100_Stnd: Growth rate as a function 
of core radius and fluid element size for Ledoux convection. 



Fig. 47. Model M25_LS_100_Stnd: Growth rate as a function 
of core radius and fluid element size for lepto-entropy fingers. 
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Fig. 45. Model M25_LS_100_Stnd: Growth rate as a function 
of core radius and fluid element size for Convection. 



and M25_LS_100Jmprv_NF are similar, respectively, to 
the corresponding graphs for models M15_LS_100_Stnd, 
M15_W_100_Stnd, and M15_LS_100Jmprv_NF. The var- 
ious modes of instability are qualitatively in the same 
regions of the core, and the growth rates of these insta- 
bilities are similar. 

We omit showing graphs of the various instability re- 
gions for the models at 200 ms as it is not clear how mean- 
ingful they are at this time. The fluid motions result- 
ing from the instabilities would have modified the back- 
ground gradients in s and Yg, thus modifying the insta- 
bilities and their growth rates. Indeed, they may already 
have done so by 100 ms. All that can be said without de- 
tailed multi-dimensional radiation-hydrodynamical sim- 
ulations is that the lepto-entropy semiconvective and 
lepto-entropy finger instabilities are likely to appear for 
for some period of time after core bounce below the neu- 
trinosphere of a proto-supernova. 
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Fig. 46. Model M25_LS_100_Stnd: Growth rate as a function 
of core radius and fluid element size for semiconvection. 



Fig. 48. Model M25_W_100_Stnd: Growth rate as a function 
of core radius and fluid element size for Convection. 



7. COMPARISONS WITH OTHER WORK 
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Fig. 49. Model M25_W_100_Stnd: Growth rate as a function 
of core radius and fluid element size for semiconvection. 



Fig. 52. — Model M25_LS_100_Imprv: Growth rate as a function 
of core radius and fluid element size for Convection. 
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Fig. 50. Model M25_W_100_Stnd: Growth rate as a function 
of core radius and fluid element size for lepto-entropy fingers. 
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Fig. 51. — Model M25_LS_100_Imprv: Growth rate as a function 
of core radius and fluid element size for lepto-entropy fingers. 



In this Section we will compare our the results of sta- 
bility analyses with the few other analyses of doubly dif- 
fusive instabilities in the literature that are relevant to 
the collapsed cores of supernova progenitors, and we will 
examine results of extant multi-dimensional radiation hy- 
drodynamical simulations of stellar core collapse to see if 
the lepto-entropy finger instability has been seen in any 
of these simulations 

7.1. Comparisons With Prior Stability Analyses 

A number of authors have recognized the poten- 
tial role of convection below the neutrinosphere in en- 
hancin g the n e utrino luminosity. Most o f this work 
^psteTrl 1197ft jBruenn et all Il979t ITivl o et alJ 

HaDT 



IColgate fc Petschekl 1198(1 ISmarr et all 
[ Lattimer fc Mazurekl [19811: [Burrows fc FrvxeTl 



Janka & Miiller 1993, 1994| iBruenn fc Mezzacappal 
1994: lKeTeTaLl [l996). was focussed on Ledoux convec- 
tion or Raleigh- Taylor overturns. ISmarr et alJ ([198 ll ) 
pointed out the possibility of doubly diffusive instabil- 
ities, in particular, neutron fingers, arising in the core 
after bounce. The Livermore group ([Wilson et al .11 198(1 
IWilson fc Mavlelll988t Il99fl have emphasized the role 
of doubly diffusive instabilities below the neutrinosphere 
and have incorporated it in their simulations. In their 
treatment of these instabilities, they have assumed that 
the rate of thermal equilibration greatly exceeds that of 
lepton equilibration, using arguments reviewed in Sec- 
tion d and on the basis of this found that material be- 
low the neutrinosphere in their core collapse simulations 
was unstable to neutron fingers in an extensive region. 
Modeling the resulting convection by a mixing length al- 
gorithm they found an enhancement of the u e luminosity 
w hich helped power ex plosions in their simulations. 

IBruenn et alJ (Il995f) a pplied the single particle anal- 
vsis of iGrossman et alJ Jl993) to stellar material with 
neutrino mediated thermal and lepton transport. On 
the basis of some numerical equilibration experiments 
they questioned the assumption that thermal equilibra- 
tion rate exceeded the lepton equilibration rate, and con- 
cluded that the neutron finger instability would proba- 
bly not be present in the collapsed cores of supernova 
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progeni tors. They also red i scover ed the the observa- 
tion of iLattimer fc Maznrekl ([198 ll) that the derivative 

( din Y e ) changes from negative to positive when the 

value of falls below a critical value, and that the core 
w as likely stable wh e n tha t occurred. 

iBruenn fc Dineval (I1996T) also applied the single par- 
ticle analysis of iGrossman et alJ ((19931 ). and extended 
the analysis by introducing the four response functions 
and the methodology for determining their values. They 
also pointed out that the rate of lepton transport sub- 
stantially exceeded the rate of thermal transport, and 
that the appearance of neutron fingers in the collapsed 
core of a sup e rnova progenitor was very unlikely. Like 
IBruenn et alJ (|1995f) . however, they applied their analy- 
sis to only a very limited region of the thermodynamic 
state space of a core and concluded that the material 
would either be semiconvective or stable, depending on 
t he value of 1/. 

iMiralles et alJ ([200(1) have examined the stability of 
the material in an evolving protoneutron star. By em- 
ploying an equilibrium diffusion approximation for neu- 
trino transport they were able to derive a dispersion 
equation (a cubic equation analogous to equation ([53[H 
involving both the growth or decay rate and the hori- 
zontal wave number or length scale of the perturbation. 
They were thus able to solve for the length scale of a per- 
turbation giving the minimum thermal and lepton gra- 
dients needed for instability or, equivalent ly, the length 
scale of a perturbation giving the maximum growth rate 
for given thermal and lepton gradients. In our work it 
was necessary to consider different length scales explic- 
itly. However, the use of an equilibrium diffusion ap- 
proximation for neutrino transport may cause some of 
the subtleties of the transport to be missed and result in 
differe nt relative values for the response functions. 

The IMiralles et alJ ((200(1 ) analysis was applied to the 
interior of a neutron star during its long term evolution. 
Their results are therefore probably not directly compa- 
rable to the results of our analysis here which pertain 
to the collapsed cores of supernova progenitors before 
the excess entropy of the shocked material has had a 
chance to be radiated away. Nevertheless we note that 
very e arly in the evolution o f the protoneutron star stud- 
ied bv IMiralles et alJ ((200(1) they find a convectively un- 
stable region near the surface surrounded by diffusively 
unstable region which they attribute to neutron fingers. 
They do not find semiconvection anywhere. This ap- 
pears to be not altogether dissimilar to our results which 
show a convectively unstable region bounded above by 
lepto-entropy fingers and below by semiconvection. The 
appearance of semiconvection below the convectively un- 
stable r egion in our models a nd not in the model exam- 
ined by IMiralles et alJ ((200(1) might be due to the fact 
that the lepton fraction in the deep interior of our mod- 
els is relatively large (so that ( j^y- ) < whereas the 

opposite might have been true in their model. Accord- 
ing to the discussion at the end of Section this could 
explai n the different results . 

The IMiralles et all ([200(1) work also included viscosity, 
the magnitude of which is uncertain and not included 
in this work, They found that viscosity, as given by 
Ivan den Horn fc van Weertl (|1981[) . did not greatly affect 



the diffusively unstable regions, but found that the con- 
vectively unstable regions to be more extended with the 
inclusion of viscosity.. 

The role of magnetic fields in the development 
of proto-neutron sta r instabilities was investigated in 
IMiralles et alJ ((20031) . Magnetic fields have a dual role: 
stabilizatio n and destabili z ation . The essential results 
presented in lMiralles et alJ ([2003D indicate that magnetic 
field strengths in excess of 10 13 G may do both: stabi- 
lize regions unstable to convection, limiting the spatial 
extent of such regions, and destabilize regions to new 
modes of instability. However, their fundamental con- 
clusion is that only very strong magnetic fields (> 10 16 
G) would significantly alter the stability properties of 
proto-neutron stars. These fields might be generated by 
dynamo action from convective activity present in the 
core prior to the escalation of field strengths. Once the 
field strength increases, convective activity is shut off. 
Thus, magnetic fields may affect both the spatial extent 
and the temporal duration of convective activity in the 
proto-neutron star. 

7.2. Comparisons With Multi- Dimensional Core 
Collapse Simulations 

Most of the extant multi-dimensional radiation- 
hydrodynamical simulations of stellar core collapse have 
employed gray neutrino transport, and it is not clear how 
much of the dynamics of the lepto-entropy finger insta- 
bility these simulations are capable of capturing. With 
this caveat in mind, we will proceed to examine the re- 
sults of these simulations that have been described in the 
literature. 

The simulations of lHerant et alJ ([1994D used a smooth 
Newtonian two-dimensional SPH (smooth particle hy- 
drodynamics) code coupled to gray neutrino transport. 
The letter transports number of energy of each neutrino 
species in a two-part scheme consisting of flux-limited 
diffusion scheme for optically thick regions and a light 
bulb approximation for optically thin regions. For their 
core collapse simulation o f a 25 M Weaver and Woosley 
IWeaver fc Wooslevl (^993) an instability below the neu- 
trinosphere developed over a time scale of 15 ms after 
bounce with maximum turnover velocities varying over 
time between 1000 and 4000 km s -1 . This unstable re- 
gion extended 15 km below the neutrinosphere which had 
a radius of 40 km at the time, with ~ 5 co nvective cells in 
a 90 degree sector. iHerant et alJ (|1994[) attributed this 
instability to the neutron finger mechanism described by 
the Livermore group. The time for the growth of the in- 
stability, the scale of the turbulence, and the region over 
which it was observed is consistent with the results dis- 
played in Figure [351 We suggest that they were seeing 
lepto-entropy fingers. 

The simulati ons of iFrverl ([19981 119991) . and 

iFrver fc Warrenl ([2002D used basically the smooth 
Newtonian dimensional SPH code couple d to gray 
neutrino transport of the lHerant et alJ (|!994f) work, but 
with improvements in the SPH code (e.g., technical 
improvements and the addition of GR effects), the 
equation of state, and the neutrino transport (the latter 
still being gray, however) . Additionally, th e simul ations 
performed perfromed by IFrver fc Warrenl ([2002D were 
three-dimensional. Unfortunately most of the discus- 
sion in these works was concerned with the entropy 
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driven convection occurring above the neutrinosphere, 
rather than con v ection below. However, Figure 2 of 
iFrver fc Warrenl ((20021) . which is a two-dimensional 
slice of the core of a 15 M© progenitor 75 ms after 
bounce, does show evidence of convective flows below 
the neutrinosphere, although it is difficult to infer 
details. 

Both the simulatio ns of Mi ller et alJ ([199 3) and 
Uanka fc Miillerl ([1996D cut out the innermost part of the 
collapsed stellar core and imposed boundary conditions, 
so fluid motions below the n eutrinosphere were not com - 
puted. In the simulations of iMezzacappa et alJ (|1998al) . 
a PPM hydrodynamics scheme was coupled to the zero 
and first moments of neutrino fields previously computed 
in spherical symmetry with a multigroup, flux-limited 
transport. While these simulations were able to cap- 
ture some convective instabilities, the imposed neutrino 
background resulted in the inability of simulations to 
compute the fluid element-to-background neutrino trans- 
port that would arise in response to differences between 
the thermodynamic state of the fluid element and the 
background. Consequently, doubly diffusive instabilities 
could not be captured by thes e simulations . 

Finally, the simulations of iBuras et alJ ((20031) used 
Newtonian PPM hydrodynamic with some general rel- 
ativistic corrections coupled to radial-ray transport via a 
multigroup, tangent ray scheme with closure provided by 
the solution of the Boltzmann equaiton on an angularly 
averaged stellar background. While details are sparse, 
the authors describe a convection that sets in below the 
neutrinosphere by 40 ms after bounce, is persistent up 
to 0.26 s after bounce when the simulations were termi- 
nated, and which slowly extended deeper into the core. 
The authors attributed this convection to Ledoux con- 
vection, but it may have been lepto-entropy fingers. 

8. CONCLUSIONS 

The fluid below the neutrinosphere of the collapsed 
core of a supernova progenitor can be driven unstable by 
thermal and lepton diffusion as well as by a gravitation- 
ally unstable stratification of the flui d. We have devel- 
oped the methodology introduced bv iBruenn Dineval 
(1996) for analyzing fluid instabilities in the presence 
of neutrino mediated thermal and lepton transfer. In 
this analysis four "response functions" are introduced to 
model the thermal and lepton equilibration of a fluid ele- 
ment perturbed in entropy and lepton fraction from that 
of the background. These response functions depend on 
the thermodynamic state of the background and the size 
of the fluid element, and are obtained by detailed neu- 
trino transport simulations. Given values of the response 
functions and the ambient gravitational field at the posi- 
tion of the fluid element, a cubic equation can be derived 
whose roots characterize the response of the fluid element 
to a perturbation in either its motion or thermodynamic 
state. From a determination of these roots, the stability 
or mode of instability of the fluid can be ascertained as 
a function of the gradients in entropy and lepton frac- 
tion. This methodology is applied both to resolve a long 
standing controversy regarding the presence of the neu- 
tron finger instability below the neutrinosphere, and to 
perform an extensive analysis of the fluid instabilities 
that do occur below the neutrinosphere for a number of 
representative post bounce supernova progenitors. 



The Livermore grou p ((Wilson et alJ Il986t 
IWilson fc Mavldll988l I1993D have argued that thermal 
equilibration proceeds much more rapidly than lepton 
equilibration in the fluid below the neutrinosphere. As a 
consequence, given the gradients in entropy and lepton 
fraction that arise in their supernova models subsequent 
to the initial collapse of the core and bounce they 
found that the neutron finger instability will arise in an 
extensive region below the neutrinosphere. Modeling by 
a mixing length algorithm the fluid motions that are 
expected to result from this instability, they found an 
enhancement of the v e luminosity that helped to produce 
successful explosions in their sophisticated spherically 
symmetric supernova simulations. Their work is very 
important as it has emphasized the potential role of fluid 
instabilities below the neutrinosphere in the supernova 
mechanism. 

We have performed cor e collapse simulations using the 
Livermore (jWilsonl 120031) equation of state. Applying 
their criteria for the neutron finger instability, we do 
find an extensive region below the neutrinosphere, on 
either side of a Ledoux unstable region, that is neutron 
finger unstable. A similar ex ercise for core collap s e sim- 
ulations computed with the IL at timer fc Swestvl ()1991t ) 
equation of state shows a similar region of neutron fin- 
ger instability, although not quite as extensive. Applying 
our methodology to the region below the neutrinosphere 
of these models we always find that lepton equilibration 
proceeds much more rapidly than thermal equilibration, 
and therefore that the neutron finger instability never 
occurs. 

On the other hand, we have also found that the cross 
response functions, viz., the tendency for a perturbation 
in entropy of the fluid element relative to the background 
to drive a net lepton flow between them, and the ten- 
dency of a lepton fraction perturbation to drive a net 
energy flow, have substantial positive magnitudes for all 
thermodynamic states, and that the ratio of the second 
to the first of these response functions is large. These 
cross response functions cause a difference in lepton frac- 
tion between the fluid element and the background to de- 
velop primarily in response to an entropy difference. This 
is in contradistinction to the case without diffusion (the 
Ledoux case) in which a lepton fraction difference devel- 
ops from a displacement of the fluid element through a 
gradient in the background lepton fraction. 

As a result of this, several new doubly diffusive instabil- 
ities appear which we have referred to as lepton-entropy 
fingers and lepton-entropy semiconvection. These insta- 
bilities are expected to present in significant regions be- 
low the neutrinospheres of post collapsed stellar cores, 
may play an important role in the supernova mechanism. 
We have analyzed these instabilities in detail. 

For each of a number of post bounce core collapse su- 
pernova models we have performed an extensive survey 
of fluid instabilities on a two-dimensional grid of fluid 
element size and radial location in the core. Our re- 
sults show a common pattern of unstable regions in these 
core collapse models. Outside of a Ledoux unstable re- 
gion is a region where the lepton fraction is small, due 
to prior shock dissociation and deleptonization, and the 
entropy gradient is positive, due to the powering up of 
the shock. The fluid is leptom-entropy finger unstable 
in this region, with maximum growth rates for fluid ele- 
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ment scales approximately 1 /20 the distance of the fluid 
element from the core center. This unstable region can 
extend up almost to the neutrinosphere, and may there- 
fore be important for the same reason adduced by the 
Livermore group for their neutron fingers. Below the 
Ledoux unstable region is region of unshocked material 
of relatively high lepton fraction. This region is unsta- 
ble to lepto-entropy semiconvection with growth rates 
favoring small scales. Reviewing the multidimensional 
core collapse simulations that have been reported in the 
literature we conclude that the lepto-entropy finger in- 
stability may have already been seen. 
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